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ABSTRACT 


Different  flow  patterns  are  encountered  in  the  calandria  of 
a  climbing  film  evaporator  when  the  feed  is  introduced  cold. 
The  lower  part  where  the  feed  is  brought  to  boiling 
temperature  is  filled  with  liquid.  The  remainder  of  the  tube 
where  evaporation  occurs  contains  liquid  and  vapor  phases. 

The  local  heat  transfer  coefficient  and  temperature 
profile  were  investigated  separately  for  the  liquid  zone 
present  in  the  lower  part  and  for  the  boiling  zone  present 
in  the  upper  part  of  a  single  tube  glass  evaporator. 

A  study  of  the  heat  transfer  in  the  liquid  zone,  where 
the  process  liquid  is  brought  to  its  boiling  temperature, 
showed  that  the  presence  of  dissolved  gases  affects 
appreciably  the  efficiency  when  the  operation  is  carried  out 
under  vacuum  as  compared  to  an  operation  at  atmospheric 
pressure . 

In  the  boiling  zone  where  two  phase-flow  is  present 
(liquid  and  vapor  phases),  a  heat  transfer  study  permits  the 
correlation  of  the  pressure  in  the  tube  with  the  evaporating 
rate  and  the  position  along  the  evaporator  tube. 

A  model  to  determine  the  position  where  boiling  begins, 
the  temperature  at  any  position  in  the  tube,  the  vapor  rate 
and  the  void  fraction  at  the  exit  of  the  evaporator  tube 
given  the  initial  conditions  (feed  rate,  feed  temperature, 
vacuum  in  the  system  and  thermal  properties  of  the  feed)  was 
set  up.  Its  validity  was  demonstrated  by  the  results 
obtained  with  sugar  solutions  and  tomato  juice.  The  model 


IV 


. 

■ 


was  also  successfully  applied  to  an  industrial  scale 


example . 
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I.  INTRODUCTION 


The  concentration  of  solutions  by  evaporation  of  part  of  the 
solvent  is  a  very  common  practice  in  the  chemical, 
pharmaceutical  and  food  processing  industries.  When  the 
material  being  handled  is  heat  sensitive,  the  problem  is 
made  more  complicated  by  the  need  to  avoid  exposure  of  the 
material  to  high  temperatures  or  a  long  contact  time  with 
the  heating  surface. 

The  use  of  a  climbing  film  evaporator  is  often  the  best 
solution,  as  high  rates  of  heat  transfer  can  be  maintained 
without  the  use  of  high  temperatures  or  long  contact  time. 

Since  1899,  when  it  was  introduced  by  Paul  Kestner,  the 
climbing  film  evaporator  has  been  used  to  a  continually 
increasing  extent  in  the  chemical  and  food  industries. 
Besides  the  short  time  of  contact  and  the  low  temperature, 
the  small  hold  up  and  relatively  simple  construction  are 
among  its  prime  advantages.  Its  use  is  not  confined  to  the 
concentration  of  aqueous  liquors,  the  same  form  of  unit  is 
also  employed  as  a  reboiler  on  fractionating  columns. 
Although  of  a  different  conception,  the  exchangers  used  for 
decompression  in  water  cooled  nuclear  reactors  show  a 
similar  flow  pattern  in  the  tubes. 

Despite  its  industrial  importance,  there  are  relatively 
few  papers  in  the  literature  in  which  precise  performance 
details  are  available,  and  the  climbing  film  evaporators  are 
designed  largely  from  experience  and  not  from  any 
established  relationship  between  the  operating  conditions, 
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the  physical  properties  of  the  liquor  and  the  consequent 
heat  flux  or  transfer  coefficient. 

This  research  was  carried  out  to  fill  this  gap  by 
establishing  relations  to  define  a  model  which  enables  one 
to  determine  the  performance  of  the  evaporator  by  control  of 
the  different  parameters  present  in  the  system. 

In  this  research  a  borosilicate  glass  climbing  film 
evaporator  was  used.  The  influence  of  the  operating  factors 
on  the  local  heat  transfer  and  temperature  profile  in  the 
liquid  zone  as  well  as  in  the  boiling  zone  is  described  and 
discussed  together  with  the  influence  of  these  same 
parameters  on  the  distance  at  which  boiling  starts  in  the 
tube  . 

The  proposed  model  is  validated  with  the  results 
obtained  with  water  and  is  successfully  tested  with  sugar 
solutions  and  tomato  juice. 

The  model  developed  for  the  single  tube  glass 
evaporator  was  tested  on  an  industrial  evaporator  made  of 
three  sections,  the  first  one  having  66  tubes,  the  second 
111  tubes  and  the  third  156  tubes  and  capable  of 
concent  rat ing  5,000  kg  per  hour  of  pineapple  juice  from 
1  3 . 5  c  B  r  i  x  to  62°Brix. 
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II.  LITERATURE  REVIEW 


A.  Climbing  film  and  vertical  tube  evaporators. 

A  review  covering  the  period  1970  to  1980  shows  that 
little  work  had  been  done  recently  in  this  area. 

The  research  on  vertical  tube  evaporators  was  generally 
conducted  in  desalination  or  wastewater  plants.  Multiple 
effect  vertical  tube  evaporators  in  desalination  plants  with 
a  capacity  of  up  to  500  tons  per  day  were  described  and 
their  performance  under  different  conditions  analysed. 

Sephton  et  al .  (1973)  and  Sephton  (1975,  1976)  showed 
that  the  addition  of  surfactant  in  sea  water  or  wastewater 
permits  a  more  stable  operation  and  an  increase  in  the  heat 
transfer  coefficient  up  to  100%  with  a  differential 
temperature  2°C  lower.  They  also  demonstrated  the 
superiority  of  upflow  vertical  tube  evaporators  versus 
downflow  vertical  tube  evaporators. 

Gel'Perin  et  a/.  (1973)  studied  the  overall  heat 

transfer  coefficient  for  boiling  mixtures  of  organic  liquids 
in  a  vertical  tube  evaporator.  They  developed  an  empirical 
equation  to  evaluate  the  overall  heat  transfer  coefficient 
in  boiling  aqueous  inorganic  salt  solutions  as  a  function  of 
the  heat  flux,  the  concentration  and  the  apparent  liquid 
level  in  a  thermosiphon  evaporator. 

Takeda  et  al .  (1972,  1973)  studied  the  effect  of  vapor 
hold  up  on  the  boiling  heat  transfer  coefficient  in  a 
natural  circulation  vertical  tube  evaporator,  but  once 
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again,  no  attempt  was  made  to  determine  the  local  heat 
transfer  coefficient  and  the  overall  heat  transfer 
coefficient  did  not  take  into  account  the  different  zones 
encountered  in  a  vertical  tube  evaporator  (liquid  not 
boiling  at  the  bottom  and  two  phase  flow  in  the  upper  part). 

Arellano  Gajardo  (1972)  used  an  apparatus  similar  to 
ours  (pyrex  glass  climbing  film  evaporator)  to  determine  the 
chemical,  nutritive  and  organoleptic  changes  in 
concen t ra t i on  of  milk,  orange  juice  and  grape  juice.  He  only 
made  a  mass  balance  study  to  determine  the  evaporative 
efficiency  of  the  climbilng  film  evaporator. 

In  previous  works  covering  the  period  1945  -  1970,  some 
authors  noticed  that  the  successive  flow  patterns 
encountered  in  the  tube  made  difficult  the  determination  of 
the  heat  transfer  coefficient  because  the  ratio  of  the 
non-boiling  zone  length  to  the  boiling  zone  length  in  the 
tube  was  a  function  of  the  feed  temperature,  feed  rate  and 
of  all  the  other  variables  in  the  system.  The  calculated 
overall  heat  transfer  coefficient  differed  from  one 
experiment  to  another.  To  avoid  this  problem,  Coulson  and 
McNelly  (1956)  introduced  the  feed  at  boiling  temperature. 
They  found  three  different  forms  of  heat  transfer  along  the 
tube  depending  on  the  feed  rate  and  the  differential 
temperature  between  the  heating  media  and  the  product.  While 
using  a  traveling  thermocouple  to  measure  the  temperature  in 
the  tube,  they  did  not  explore  the  full  capacity  of  their 
apparatus  and  kept  calculating  overall  heat  transfer 
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coefficients  based  on  the  average  temperature  in  the  tube. 

Guerrieri  and  Talty  (1956)  installed  thermocouples 
every  6  inches  in  a  6  feet  long  vertical  tube  evaporator  and 
studied  the  local  heat  transfer  to  organic  liquids,  but 
unfortunately  they  developed  their  study  from  the  Lockhart 
and  Martinelli  (1949)  pressure  drop  correlation  for 
two-phase  flow.  This  correlation  was  developed  for 
isothermal  two-phase,  two-component  flow  in  pipes  without 
phase  change,  which  does  not  correspond  to  the  experimental 
conditions  encountered  by  Guerrieri  and  Talty. 

B.  Two- phase  flow 

A  literature  review  on  two-phase  flow  shows  that  the 
prediction  of  heat  and  mass  transfer  when  a  liquid  phase  and 
a  gaseous  phase  are  present  in  a  tube  is  difficult  because 
of  the  number  of  parameters  involved.  The  present  knowledge 
of  heat  transfer  in  such  flow  systems  is  limited  and  results 
obtained  for  one  set  of  conditions  are  of  little  value  in 
predicting  performance  at  a  different  set  of  conditions. 

Lunde  (1961)  tried  to  find  a  correlation  for  two-phase, 
two-component  heat  transfer  on  the  basis  of  a  physical 
model.  He  took  the  data  obtained  by  different  authors  and 
correlated  them  for  the  supposed  different  flow  patterns 
encountered  (stratified  flow,  slug  flow,  annular  flow);  in 
most  cases  the  deviation  of  the  experimental  data  from  the 
proposed  correlation  is  up  to  30%.  Sometimes  points  are 
outside  the  30%  deviation  zone.  This  model  permits  only  the 
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evaluation  of  an  overall  heat  transfer  coefficient.  This 
model  is  probably  acceptable  when  the  change  of  phase  in  the 
exchanger  is  kept  at  a  minimum  to  ensure  that  the  ratio  of 
gas  to  liquid  does  not  vary  too  much  and  that  the  assumed 
flow  pattern  stays  the  same  all  along  the  tube,  which  is  not 
the  case  in  our  study. 

Lunde  (1961)  showed  in  his  study  that  the  critical 
Reynolds  Number  for  the  transition  between  viscous  and 
turbulent  regions  is  smaller  for  two-phase  flow  than  for 
one-phase  flow.  This  phenomenon  was  confirmed  by  Kao  et  al . 
(1972)  in  their  analysis  of  the  stability  of  film  boiling  on 
vertical  surfaces. 

A  complete  analysis  of  one-dimensional  two-phase  flow 
was  made  by  Wallis  (1969).  A  review  of  all  the  possible 
problems  encountered  in  one-dimensional  two-phase  flow  was 
given  and  equations  proposed  by  different  authors  were 
discussed.  When  no  phase  change  occurs  in  the  tube,  the 
validity  of  the  proposed  equations  which  are  based  on 
continuity  and  momentum  is  without  equivocation;  but  when  it 
comes  to  flow  with  phase  change,  the  different  authors  do 
not  agree  on  the  repartition  of  the  force  for  momentum 
increase  due  to  phase  change  per  unit  volume  between  the 
liquid  phase  and  the  vapor  phase.  This  repartition  appears 
to  be  function  of  the  process,  the  hydrodynamic  mechanism 
and  the  rate  of  vaporisation. 

Malnes  (1975),  in  his  analysis  of  critical  two-phase 
flow  based  on  non-equilibrium  effects,  showed  that  in  real 
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systems,  two-dimensional  effects  occur  and  affect  the  values 
predicted  by  a  one-dimensional  theory  by  an  increase  of  the 
exit  pressure  and  a  decrease  of  exit  velocities.  To  obtain 
results  in  accordance  with  reality,  Malnes  (1975)  had  to 
introduce  a  correlation  for  the  slip  ratio  between  the  vapor 
phase  and  the  liquid  phase  and  a  correlation  for  flashing 
occur i ng • along  the  tube  of  the  evaporator. 

c.  Climbing  film  evaporators  in  the  food  industry. 

In  food  processing,  evaporation  implies  a  unit 
operation  where  the  concentration  of  a  solution  is  increased 
by  driving  off  a  portion  of  the  solvent  in  the  form  of 
vapor.  In  most  cases  the  solvent  vaporised  is  water  and  the 
concentrated  solution  is  the  desired  product. 

The  problem  is  far  more  complex  than  the  evaporation  of 
most  chemical  products,  as  the  end  product  is  not  only 
subject  to  physical  and  biological  testing  but  to 
organoleptic  analysis.  Such  hard-t o-def ine  parameters  as 
taste,  odor,  color  and  texture  must  be  of  prime  concern  in 
the  process  design  as  well  as  in  the  equipment  design  and 
operat i on . 

It  is  essential  that  the  evaporation  of  water  from  the 
product  by  heat  shall  be  mild  enough  to  preserve  the 
nutritive  value  of  the  product  and  prevent  browning 
reactions,  either  by  caramelisation  of  sugars  or  by  reaction 
between  protein  and  sugars.  Also,  undue  aenaturation  of  the 
proteins  and  destruction  of  vitamins  by  excess  or  prolonged 
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heating  is  undesirable. 

The  oldest  of  the  evaporators  is  the  shell  and  tube 
type.  The  first  real  advancement  in  evaporating  equipment 
was  made  in  1813  when  E.  C.  Howard  invented  the  vacuum  pan 
which  reduced  the  boiling  temperature  by  about  one-half.  As 
stated  earlier,  the  next  major  step  occurred  in  1899  when 
P.  Kestner  invented  the  vertical  tube  climbing  film 
evaporator . 

Unf ortunatly ,  the  climbing  film  evaporator  as  a  unit 
has  not  found  much  favour  in  the  food  industry.  Many 
evaporators,  if  not  classified  as  true  climbing  film, 
contain  a  section  where  climbing  film  evaporation  is 
encountered.  It  is  useful  to  review  these  evaporators  and 
their  use  in  the  food  industry. 

1.  Climbing  film  evaporators. 

Sometimes  it  is  necessary  to  produce  from  a  very  weak 
liquor  a  highly  concentrated  finish  product.  A  single  pass 
climbing  film  evaporator  will  not  solve  this  problem.  One 
solution  is  to  recirculate  the  liquor  through  the 
evaporator.  Another  solution  is  to  pass  the  liquor  in  series 
through  two  or  more  evaporator  units.  This  is  known  as  the 
mult i -c i rculat ion  climbing  film  evaporator  (Slade,  1967). 

It  includes  two  or  more  calandrias.  These  calandrias 
are  formed  by  dividing  a  simple  steam  shell  into  two  or  more 
sections  of  tubes,  each  section  having  fewer  tubes  than  the 
preceding  one  and  being  provided  with  a  separator  at  the 
top.  The  flow  of  liquor  passing  through  each  section  and 
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each  separator  in  series  progressively  diminishes  in 
quantity,  but  since  the  number  of  tubes  is  progressively 
reduced,  an  adequate  flow  is  provided  in  each  section  to 
maintain  a  climbing  film  condition.  Such  an  evaporator  has 
been  used  for  handling  heat  sensitive  liquors  such  as 
coffee,  fruit  juices,  milk,  tea,  whey  and  wine. 

The  calandria  short  vertical  tube  evaporator  as 
described  by  Moore  and  Hesler  (1963)  has  been  used  for  over 
a  century  in  the  sugar  industry  but  cannot  be  considered  as 
a  true  climbing  film  evaporator,  the  tubes  being  too  short 
to  obtain  a  true  film  boiling  zone. 

The  sloping  calandria  evaporator  described  by  Brennan 
et  al  .  (  1  969)  looks  like  a  recirculating  climbing  film 
evaporator  tilted  in  an  inclined  position.  It  was  widely 
used  in  the  dairy  industry  up  to  the  mid-fifties  when  it  was 
progressively  replaced  by  more  sophisticated  evaporators 
(Scott ,  1964). 

Liquid  food-stuffs  very  sensitive  to  heat  can  be 
evaporated  at  very  low  temperature  in  a  low  temperature 
evaporator  (Slade,  1967).  Such  an  evaporator  is  of  the 
recirculating  climbing  film  type  and  is  used  under  high 
vacuum.  The  boiling  point  of  the  liquor  being  of  the  order 
of  15°C,  the  temperature  of  the  condenser  water  usually 
available  for  other  types  of  evaporators  is  too  high  to 
condense  the  vapor  formed  at  such  a  low  pressure;  this  is 
why  an  auxiliary  fluid  (ammonia)  is  used.  It  works  on  the 
principle  of  the  heat  pump;  the  heating  fluid  is  vaporized 
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to  cool  the  condenser,  then  recompressed  through  a 
compressor  to  supply  heat  to  the  calandria. 

The  low  temperature  evaporator  was  widely  used  in  the 
concentration  of  orange  juice  following  the  Second  World  War 
but  was  replaced  by  single  pass  evaporators,  where  the 
advantage  of  the  low  temperature  was  overcome  by  the 
decrease  in  residence  time  leading  to  similar  or  even  better 
organoleptic  qualities  of  the  concentrated  juice. 

2.  Climbing  and  falling  film  evaporators. 

A  combination  of  the  climbing  and  falling  film  may  be 
used  where  a  high  percentage  of  evaporation  is  required  and 
when  the  concentrated  liquor  tends  to  be  viscous;  as  the 
most  concentrated  liquor  is  formed  in  the  falling  film 
section,  the  flow  is  improved  by  the  gravity  force.  It  is 
generally  used  as  a  finisher  evaporator  following  a  climbing 
film  evaporator  where  the  greatest  part  of  the  evaporation 
has  been  effected.  It  has  been  used  to  produce  orange  and 
tomato  juices. 

3.  Plate  evaporators. 

The  plate  evaporator  is  well  described  by  Carter  and 
Kraybill  (1966).  It  differs  from  the  conventional  evaporator 
in  that  the  tubular  calandria  is  replaced  by  a  series  of 
plates.  The  plates  are  arranged  so  that  the  liquid  passes 
through  a  rising  film  section,  then  down  through  a  falling 
film  section  with  each  product  plate  sandwiched  between  two 
steam  plates.  A  very  high  surface  to  volume  ratio  is 
achieved,  thereby  reducing  the  holdup  time. 
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Sanitary  type  stainless-steel  centrifugal  pumps  are 
generally  employed  for  liquor  feed,  and  a  velocity  up  to 
80  m.s'1  can  be  achieved  in  the  last  plates  of  the 
evaporator  ( Papadopoulos  and  Samuel,  1965). 

Among  the  other  advantages  of  the  plate  evaporator  are 
its  compactness,  the  possibility  of  increasing  the  capacity 
rapidly  by  adding  plates  and  the  good  access  for  cleaning  by 
quick  di sconnect ion  of  the  plates.  A  further  obvious 
advantage  is  that  should  any  plate  failure  occur,  a 
replacement  can  be  made  immediately  with  a  minimum  loss  of 
production  time. 

Lawler  (1960)  and  Papadopoulos  and  Samuel  (1965) 
discuss  the  use  of  this  evaporator  in  the  concentration  of 
orange  juice  and  beer  wort.  It  is  widely  in  use  all  over  the 
world  to  concentrate  a  great  variety  of  products  such  as 
skim  milk,  whole  milk,  buttermilk,  whey,  ice  cream  mix, 
apple  juice,  orange  juice,  pear  juice,  lemon  juices,  apple 
pectin,  chicken  broth,  gelatin,  jam  and  table  jelly. 

4.  Expanding  flow  evaporators. 

This  type  of  evaporator  is  almost  exclusively  used  in 
the  dairy  industry.  It  is  well  described  by 
Brennan  et  al .  (1969).  In  this  device  milk  and  steam  flow 

through  alternate  passages  in  a  manner  similar  to  that  in  a 
plate  evaporator.  The  plates  are  replaced  by  thin,  inverted 
stainless  steel  cones,  gasketed  to  prevent  leakage.  Milk 
enters  through  a  central  spindle  at  the  base  of  the  cone 
assembly  and  is  fed,  via  nozzles,  to  the  heated  cone  spaces. 
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Flowing  upwards  and  outwards  over  the  steam-heated  surfaces, 
under  vacuum  conditions  the  milk  rapidly  attains  its  boiling 
point.  A  high  velocity  vapour-liquid  stream  leaves  the  cone 
assembly  tangentially,  ensuring  a  good  separation. 


III.  EXPERIMENTAL 


A.  Equipment. 

The  evaporator  is  a  Q.V.F.  CFE1  climbing  film 
evaporator  made  of  borosilicate  glass.  A  schematic  diagram 
of  the  evaporator  is  presented  in  Figure  1. 

It  is  composed  of  three  parts: 

1.  The  evaporator  part  which  is  made  of  a  calandria  tube  3m 
long,  3cm  inside  diameter  and  1.53mm  wall  thickness. 

This  calandria  is  enclosed  in  a  steam  jacket  made  of  a 
2.70m  long,  5cm  bore  glass  pipe  which  can  withstand 
steam  pressure  of  3  bars  absolute.  It  has  three  branches 
for  steam  inlet,  condensed  steam  and  vent  connections. 

2.  The  inside  evaporator  tube  leads  from  the  calandria  to 
the  condenser  via  a  cyclone  separator  which  separates 
the  entrained  liquid  from  the  vapor.  The  liquid  outlet 
of  the  separator  is  connected  directly  to  the 
concentrate  receiver.  The  separator  and  the  concentrate 
receiver  form  the  second  part  of  the  apparatus. 

3.  The  third  part  of  the  evaporator  is  composed  of  the 
vapor  condenser  directly  connected  to  the  vapor  exit  of 
the  cyclone  separator.  The  bottom  part  of  the  condenser 
leads  to  two  5  liters  condensate  receivers. 

Each  condensate  receiver  is  connected  to  the  vacuum 
line.  The  two  receivers  are  separated  by  a  valve  to  enable 
the  condensate  to  be  removed  during  the  operation  of  the 
evaporator.  When  the  lower  one  is  full  of  condensed  water, 
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Pressure  tap 


Fig.  1.  Schematic  diagram  of  the  evaporator. 


15 


it  can  be  isolated,  vented  to  atmospheric  pressure,  emptied 
and  vacuum  can  be  reapplied  without  interruption  of  the 
experiment . 

The  unit  is  connected  to  a  Nash  Hytor  vacuum  pump 
model  673  capable  of  maintaining  a  pressure  of  160mm  of 
mercury  in  the  system.  The  vacuum  can  be  adjusted  for 
operation  and  released  at  the  end  of  the  operation  by  means 
of  a  vent  cock  on  the  main  line. 

Seven  thermocouples,  type  T  (copper  constantan),  are 
placed  along  the  axis  of  the  calandria  tube.  The  wires  are 
0.46mm  in  diameter  and  a  common  constantan  wire  is  used  to 
reduce  the  number  of  wires  in  the  tube. 

Thermocouple  1  is  located  under  the  steam  jacket 
connection  and  records  the  feed  temperature.  Thermocouples 
2,  3,  4,  5  and  6  are  located  at  0.15m,  0.45m,  0.95m,  1.45m 

and  1.95m  from  the  condensed  steam  vent  which  is  taken  as 
origin  to  measure  the  position  along  the  heat  exchanger 
calandria.  Thermocouple  7  is  located  just  above  the  steam 
jacket  and  records  the  exit  temperature. 

Other  thermocouples  are  located  on  the  top  and  the 
bottom  of  the  vapor  condenser,  at  the  inlet  and  at  the 
outlet  of  the  cooling  water  line  of  the  condenser.  Another 
thermocouple  is  kept  in  a  bath  of  ice  in  water  and  is  used 
as  reference  to  zero  the  recorder. 

The  twelve  thermocouples  are  connected  to  a  twelve 
point  multipoint  recorder  (Honeywell  Electronik  model  112) 
calibrated  in  degrees  centigrade  with  a  range  0  to  1  00  0  C . 
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A  pressure  tap  located  between  the  cyclone  separator 
and  the  vapor  condenser  is  connected  to  a  mercury  manometer 
to  obtain  an  accurate  reading  of  the  vacuum  in  the  system. 

The  steam  is  supplied  through  a  variable  pressure 
control  valve  which  adjusts  the  steam  pressure  to  the 
desired  value.  The  feed  rate,  when  using  a  clear  solution, 
is  measured  through  a  set  of  two  rotameters  with  a  maximum 
capacity  of  3kg/min  of  water  at  20°C. 

B .  Method . 

In  normal  use  of  the  evaporator,  different  zones  are 
present  in  the  calandria  tube  (liquid  zone  at  the  lower  part 
and  film  boiling  zone  at  the  upper  part).  Because  in  each 
zone  the  heat  transfer  is  controlled  by  a  different 
mechanism,  the  study  was  divided  into  three  parts: 

1.  Experimentation  with  the  liquid  phase  filling  completely 
the  tube  with  no  evaporation.  The  evaporator  is  used  as 
a  single  tube  heat  exchanger.  This  gives  the  correlation 
of  the  local  heat  transfer  in  the  liquid  zone. 

2.  Experimentation  with  the  liquid  entering  the  evaporator 
at  its  boiling  temperature  to  obtain  a  two-phase  flow 
pattern  all  along  the  tube  and  develop  a  model  for  the 
climbing  film  ebullition  zone. 

3.  In  the  third  part,  a  model  is  proposed  which  by 
combination  of  the  relationships  obtained  in  the  two 
previous  sections  permits  the  prediction  of  the 
performance  of  the  evaporator  in  normal  use  (distance  at 
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which  boiling  starts,  vapor  rate  at  the  exit  and 
temperature  profile  along  the  tube)  knowing  the  inlet 
condi t i ons . 

The  correlations  governing  the  liquid  zone  were 
determined  from  data  obtained  using  water  as  a  medium.  In 
this  section,  the  evaporator  was  run  at  atmospheric  pressure 
and  the  water  was  fed  by  means  of  a  variable  flow  positive 
pressure  pump.  The  lower  flow  rate  is  limited  by  the  water 
reaching  boiling  temperature  at  the  exit  of  the  evaporator 
tube  while  the  maximum  rate  is  limited  by  the  capacity  of 
the  rotameter. 

Correlations  in  the  film  boiling  zone  were  determined 
using  water  as  medium.  The  evaporator  was  run  under 
controlled  vacuum  and  the  feed  temperature  was  adjusted  to 
the  boiling  temperature  corresponding  to  the  pressure  in  the 
system  by  means  of  an  heat  exchanger  mounted  prior  to  the 
evaporator  inlet. 

The  feed  is  introduced  under  a  pressure  gradient 
existing  between  the  container,  where  it  is  stored  at 
atmospheric  pressure,  and  the  vacuum  in  the  system.  The  rate 
is  controlled  by  means  of  a  needle  valve  placed  on  the  feed 
line. 

The  validity  of  the  model  developed  in  the  third  part 
was  tested  using  water,  sucrose  solutions  of  different 
concentration  and  commercial  tomato  juice. 

A  typical  run  with  the  evaporator  is  as  follows: 

1.  The  steam  pressure,  vacuum  and  feed  rate  are  set  at  the 
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desired  values. 

2.  After  a  few  minutes  when  a  steady  state  has  been  reached 
(the  temperature  profile  along  the  tube  does  not  change 
any  more),  the  level  of  concentrate  in  the  concentrate 
receiver  is  noted  and  the  time  is  set  at  zero. 

3.  When  either  the  concentrate  receiver  or  the  condensate 
receiver  are  almost  full,  the  level  in  the  concentrate 
receiver  is  recorded  again  at  time  t.  The  concentrate 
receiver  having  been  calibrated,  the  rate  of  concentrate 
collection  is  evaluated. 

4.  The  rate  of  vapor  formation  at  the  exit  of  the 
evaporator  is  obtained  by  difference  between  the  feed 
rate  and  the  concentrate  rate. 
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IV.  THEORY 


A.  Liquid  zone. 

When  there  is  only  the  liquid  phase  present  in  the 
tube,  considering  a  tube  section  of  length  dz ,  the  heat  flow 
from  the  steam  to  the  liquid  can  be  evaluated  as  follows 
(refer  to  Figure  2 ) : 

q  =  W.Cp.dT  =  Uz .dA. (Ts  -  T)  (1) 

with  dA  =  7T.D.dz  (2) 

Solving  for  Uz  we  obtain: 


W.Cp  dT 

Uz.dz  =  -  -  (3) 

7T.D  (Ts  ~  T) 


D imen s i onal i za t i on  of  the  equations: 


z  dz 

Z*  =  -  dZ*=  —  (4) 

L  L 


(Ts  -  Ti  )  (Ts  -  Tj  ) 

T*  =  -  dT*  =  -  dT  (5) 

(Ts  -  T)  (Ts  -  T)  2 
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Condensate 


2.  Cross  section  of  the  evaporator  tube  in  the  liquid 
zone . 


Fig 


Equation  3  becomes: 


W.Cp  dT* 

Uz  .  dZ  *  =  -  -  (6) 

7T.D.L  T* 


Because  the  specific  heat  of  liquid  food  varies  by  less  than 
1%  over  the  temperature  range  0°C  to  100°C,  we  assume  that 
W .  Cp/ 7T.  D  .  L  is  constant. 

Integration  of  Equation  6  leads  to: 


z  * 

W.Cp  T* 

Uz .dZ* 

In 

Z  o  * 

7T.D.L  T0* 

(Ts  "  Ti  ) 

with  T0*  =  -  =  1 

(Ts  -  T-t  ) 


Z*  W.Cp 

Uz  .dZ*  =  -  In (T* ) 

Z  o  *  7T.D.L 


(7) 


For  a  given  flow  rate  and  a  known  initial  temperature, 
ln(T*)  is  only  a  function  of  the  position  Z*  in  the  tube. 
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Equation  6  can  be  rewritten  as: 


W.Cp  dln(T*) 

Uz  = -  (8) 

7T.  D .  L  dZ* 


Let  us  introduce  the  Peclet  Number  which  is  the  product  of 
the  Reynolds  Number  by  the  Prandlt  Number. 


4.W  /x.cp 

Re  =  -  and  Pr  =  - 

7T.jU.D  k 


4 .W.Cp 

Pe  =  Re  .  Pr  =  -  (9) 

7T.D.  k 


The  initial  Peclet  Number  (at  the  inlet  of  the  evaporator 
tube)  is  defined  as  follows: 

4 .W.Cp 

Pe,  =  -  (10) 

7T .  D  .  k  i 


Equation  8  becomes: 


Pei  dln(T*) 

Uz  =  —  .  ki  - 

dZ* 


4  .L 


(11) 
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Rearrangement  of  Equation  11  leads  to: 


U  z . L  Pei  din ( T* ) 

kj  4  dZ* 

Uz .L/k;  can  be  defined  as  a  local  Nusselt  Number  (Nuz) 
characteristic  of  our  evaporator  length  and  overall  heat 
transfer  coefficient  at  position  z. 

Equation  11  becomes: 


Pei  dln(T*) 

Nu  z  =  -  -  (12) 

4  dZ* 

With  in(T*)  only  a  function  of  Z*  and  Pe j .  Therefore  the 
heat  transfer  in  the  liquid  zone  should  follow  an  equation 
of  the  form: 


Nuz  =  f (Pei ,  Z*) 


(13) 


This  relation  is  to  be  determined  empirically  from 
experimental  data.  This  will  be  done  in  the  Results  Section. 
When  Equation  13  is  known,  it  will  be  used  to  predict  the 
temperature  in  the  liquid  zone  using  Equation  7. 
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B .  Boil ing  zone. 

1.  Pressure  along  the  tube. 

We  assume  that  the  vapor  phase  is  always  saturated. 
Using  the  equations  proposed  by  Malnes  (1975)  we  define 
the  continuity  and  momentum  equations  as  follows: 

Continuity: 


dG  i 

-  =  -  \p  (14) 

dz 


dG2 

-  =  \l/  (15) 

dz 


Momentum : 


dP  d 

0  =  —  -  —  ( G , . V ,  +  G2.V2)  -  F  -  gr  (16) 

dz  dz 


where 

F  =  Friction  term. 


gr  =  Gravity  head. 


25 


We  define  the  slip  ratio  as  follow: 


V2 

S  =  —  (17) 

V, 


Therefore  we  have,  assuming  that  S  stays  constant  (refer  to 
the  Discussion  of  Results  Section): 


Qj 

< 

to 

dV, 

=  s  — 

(18) 

dz 

dz 

Levy  (1960)  suggested  that,  if  evaporation  is  rapid  enough 
(which  is  the  case  in  our  study),  inertia  and  phase  change 
terms  are  usually  much  larger  than  friction  and  body  terms. 
So  we  assume  F  =  0  and  Equation  16  becomes: 


dP 

dG  , 

dV, 

dG  2 

dV2 

—  +  V, 

+  G, 

+  v2 

+  G  2  +  gr  -  0 

(19) 

dz 

dz 

dz 

dz 

dz 

The  mass  fluxes  are  defined  as  follows: 


G,  =  p  ,  .  V ,  (  1  -  a)  (20) 


G  2  Q  2  •  ^  2  •  ^ 


(21  ) 
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Development  of  Equations  14  to  21  leads  to: 


dP  dV,  dV, 

—  -i/'.V,  +p,  .V,  (  1  -  a)—  +S.V,  .i//+p2  .S2  .V,  .a—  +  gr  =  0  (22) 

dz  dz  dz 


Differentiation  of  Equations  20  and  21  with  respect  to  z 
gives,  respectively: 


dG, 

dV, 

dp , 

da 

-  -  \ 

ii 

i 

R 

• 

"tT 

+  V,  ) 

~  P  i  •  v , 

(23) 

dz 

dz 

dz 

dz 

and 

dG  2 

dV, 

dp  2 

da 

-  ; 

j/  -  S  .  a  .  ( p  2 

+  V,  ) 

+  S  .p  2  .  V ,  — 

(24) 

dz 

dz 

dz 

dz 

Elimination  of  da/dz  between  Equations  23  and  24  and  solving 
for  dV,/dz  gives: 


dV,  (p,  -  S.p2) 

-  =  ^  - 


a  dp  2 

V,  ( - 

p  2  dz 


(  1  -  a  )  dp  , 

- ) 

P. 


dz 


S  •  p^  i  •  P  2 


+ 


dz 


(25) 
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Equations  22  and  25  give: 


dP  (p ,  -  S ,p2  )  a  dp  2 

—  +V,.^.  (p  -  +  S  -  1  )  -p  .  V ,  2  ( - 

dz  S .p, ,p,  p  j dz 


(  1  -  a  )  dp  , 
p  i  dz 


) +  gr  =  0 


(26) 


with  p  =  a .p  2 . S  2  +  (1  -  a )p , 


(27) 


dp,/dz  and  dp2/dz  are  as  follows: 


dp , 

dp , 

b  T 

dp , 

dP 

dz 

dT 

dz 

dP 

dz 

dp  2 

dp  2 

dP 

dp  2 

dT 

dz 

dP 

dz 

dT 

dz 

Because  the  liquid  phase  is  incompressible  we  have: 
dp,/dP=0.  The  temperature  does  not  vary  too  much  in  the 
boiling  zone  and  dp,/dT.dT/dz  can  be  neglected,  so  we  assume 
dp , /dz=0 . 

The  temperature  of  saturated  vapor  is  only  a  function 
of  the  pressure;  therefore  p2  is  only  a  function  of  the 
pressure . 


We  have  now: 
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dp, 

-  =0  (28) 

dz 


dp 2  bp2  dP 
dz  dP  dz 


(29) 


Introduction  of  Equation  29  in  Equation  26  gives: 


dP  a  dp2  (p,  “  S.p2) 

—  (1  -  p.V,2  —  - )  =  -  V,  -  +  S-1)  -  gr  (30) 

dz  p  2  dP  S .p , .p  2 

Along  the  saturation  curve  of  water  the  absolute  temperature 
and  the  pressure  can  be  correlated  by: 


T  =  c  .  P  n 

T  (K) 

P  (Pa) 


with 


/  c  =  166.766 


<  and 


V  n  =  0.06978 


(31  ) 


For  ideal  gases  we  have: 


P  =  p2.R.  T 


(32) 


where  R  is  the  gas  constant  defined  as  follows: 
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R  =  R/M  where  M  is  the  molecular  weight  and  R  is  the 
universal  gas  constant  and  is  equal  to 
8314.5  J . kmole ~ 1 . K" 1 

then  P  =  p2.R.c.  Pn  —  p2  =  (  1/7?.c  )  .  P1 

bp  2  1  -  n 


P- 


(33) 


5p  R .  c 


The  gravity  head  is  defined  as  follows: 


gr  =  g . (a.p2  +  p, ( 1  -  a)) 


(34) 


Equation  30  becomes: 


(p,  -  S  .p2  ) 


V, 


+  S  ~  1)  ~  g .  ( a  .  p  2  +p,(l  -  a)) 


dP 


S  .  P  i  .  P  2 


(35) 


az 


a  (1  -  n ) 


1  -  p.V,  2 


p  2  R.c 


Let  us  introduce  the  dimensionless  length  Z*  and  the  rate  of 
evaporation  of  water  per  unit  volume: 
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dW2 

- 

7T.  H  2  . L . dZ* 


(36) 


Introduction  of  Equation  36  in  Equation  35  gives: 


V,  (p,  -  S.p2)  dW 2 

-  - ( p  -  +  S  -  1  )  -  -  g  .  ( a  .  p  2  +  p ,  (  1  -  a  )  )  L 

dP  7T.R2  S.p,  .p2  dZ* 

-  =  -  (37) 

dZ*  a  ( 1  -  n ) 

1  -  p  .  V ,  2  —  -  P  - n 

p2  R.c 

The  pressure  can  be  obtained  in  the  boiling  zone  by  solving 
Equation  37  if  dW2/dZ*  and  a  are  known. 

2.  Void  fraction. 

a  is  not  known  but  it  can  be  estimated  with  the 
modified  Darmouth  correlation  for  vertical  upward  cocurrent 
annular  flow  as  proposed  by  Turner  (1966): 


j  2  *  j  1  * 

- - =  !  (38) 

1  -  3  .  1  (  1  -  a )  3 .  1  (  1  -  a ) 


Where  j,*  and  j2*  are  the  dimensionless  volumetric  fluxes 
defined  below. 


j  i  Vp”- 

j  i  *  = 


j  2  -\J  P  2 

j2*  = 


V g .D(p,  -  p2 ) 


and 


Vg  .D(p,  -  p2  ) 
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3.  Heat  transfer  coefficient. 

The  heat  balance  at  position  z  is  defined  as  follows 
( refer  to  Figure  3 ) : 


q  =  U  z  .  7T.  D .  dz  (Ts  -  T)  =  X.dW2  +  (W,.Cp,  +  W2.Cp2).dT  (39) 

T  is  the  temperature  of  the  liquid  at  the  interphase.  It  is 
slightly  higher  than  the  vapor  temperature  due  to  the  effect 
of  the  boiling  point  rise  for  solutions  other  than  pure 
water . 

Equation  39  leads  to: 


dW2  dT 

Uz.7T.D(Ts  -  T)  =X -  +  (W,  .Cp,  +  W2.Cp2)  —  (40) 

dz  dz 


Introduction  of  the  dimensionless  length  gives: 


dW2  dT  dP 

U  z  .  7T.  D .  L  ( Ts  -  T)  =  X -  +  (W,  .Cp,  +  W2.Cp2)  —  —  (41) 

dZ*  dP  dZ* 


Recall  Equation  31: 


Ts  = 


c  .Psn 


and 


T 


c  .Pn 
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Condensate 


Fig.  3.  Cross  section  of  the  evaporator  tube  in  the  boiling 
zone . 
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By  differentiation  with  respect  to  P  we  obtain: 


dT 

—  =  n .c  .  P  n~ 1  (42) 

dP 


Introduction  of  Equation  42  in  Equation  41  leads  to: 


U: 


7T.D.L(Psn  -  Pn) 


X  dW2  dP 

-  -  +  (W,.Cp,  +  W2 .Cp2 )n.P"' ’ — 

c  dZ*  dZ* 


(43) 


Equation  43  will  be  solved  numerically  together  with 
Equation  37.  This  will  give  values  of  the  local  heat 
transfer  coefficient  for  the  boiling  zone.  These  values  of 
Uz  ,  after  being  correlated,  will  be  used  in  the  testing  of 
the  model  to  evaluate  dw2/dz*. 


■ 
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V.  RESULTS 


a.  Liquid  zone. 

1.  Experimentation  under  atmospheric  pressure. 

Seven  different  flow  rates  of  water  were  used,  the 
system  was  kept  at  atmospheric  pressure  and  the  steam 
pressure  set  at  1.357  10s  Pa  absolute.  When  a  steady  state 
was  reached,  the  temperature  was  recorded  for  a  period  of  at 
least  15  min.  An  average  value  of  the  temperature  read  by 
each  of  the  7  thermocouples  over  this  period  was  taken.  The 
temperature  fluctuation,  which  was  nil  for  thermocouples  1 
and  2,  increased  with  the  distance  to  reach  a  value  of  1°C 
at  thermocouple  7  (exit  temperature) 

The  temperature  profile  along  the  tube  as  a  function  of 
the  distance  z  is  shown  for  different  flow  rates  in 
Figure  4.  From  Figure  4  it  appears  that  the  lower  part  of 
the  exchanger  tube  does  not  seem  to  be  effective;  this  was 
probably  due  to  the  fact  that  the  flow  was  laminar  prior  to 
the  heat  transfer  zone  and  the  the rmocouples  placed  at  the 
center  of  the  tube  did  not  record  any  temperature  change. 
After  a  short  length  of  heating,  the  tube  being  vertical, 
convection  occurred  and  the  velocity  profile  was  no  longer 
parabol ic . 
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4.  Temperature  profile  in  the  liquid  zone. 
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Recall  Equation  12: 


Pe i  dln(T*) 

Nu:  =  -  - 

4  dZ* 


Looking  at  Figure  5,  showing  ln(T*)  versus  Z*,  we  can  see 
that  for  a  given  flow  rate  ln(T*)  can  be  approximated  by  a 
second  degree  polynomial  in  Z* : 


ln(T*)  =  A  +  B . Z  *  +  C.Z*2 


(44) 


The  coefficients  A,  B  and  C  for  each  flow  rate  are  given  in 
Table  1.  Z0*  is  the  root  of  Equation  44  when  T*=1  in  the 
interval  [0,  1],  It  gives  the  distance  at  which  the  heat 

transfer  apparently  starts  and  reflects  entrance  effects. 

The  inlet  conditions  are  known,  therefore  the  initial 
Peclet  Number  (Pe-,  )  is  known. 

Figure  6  shows  the  values  of  A,  B  and  C  versus  1/Pei . 
The  lines  represent  the  best  least  squares  fit  of  A,  B  and  C 
function  of  1/Pei  for  a  regression  of  degree  1  for  A  and 
degree  2  for  B  and  C.  The  squares  of  multiple  correlation 
coefficients  are,  respectively,  0.992  for  A,  0.9995  for  B 


and  0 . 9985  for  C . 


- 

r'.rk^t  4 
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Fig .  5 . 


Liquid  zone:  In  of  the  dimensionless  temperature 
T*  versus  the  dimensionless  position  Z* 


kg  s  1 O2 
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Coefficients  A,  B  and  C  from  Equation  44  vs  1/Pei  . 
The  curves  represent  the  best  least  squares 
fit. 


Fig.  6 . 
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A  =  0.01698  -  494.9/Pei  (45) 


B  =  -0.458  +  9445/Pei  +  5.22  106/Pei2  (46) 


C  =  0.238  +  755/Pei  -  135  1 0 1 /Pe i  2 


(47) 


Now  Equation  12  becomes: 


Pe  i 

Nuz  =  - (B  +  2C . Z  * )  (48) 

4 


This  equation  permits  the  evaluation  of  the  local  heat 
transfer  coefficient  in  the  liquid  zone  when  no  gas  bubbles 
are  present. 

2.  Experimentation  under  vacuum. 

When  using  the  evaporator  under  vacuum  with  water  from 
the  tap,  we  found  that  the  temperature  increases  faster  and 
earlier  than  predicted  under  atmospheric  pressure. 

Looking  at  the  liquid  zone,  we  can  see  a  multitude  of 
tiny  bubbles  of  gas  which  are  released  from  the  water  under 
vacuum.  If  we  degassed  the  water  before  evaporation,  the 
temperature  profile  under  vacuum  remained  the  same  as  for 
atmospheric  pressure,  as  shown  in  Figure  7.  The  heat 
transfer  in  the  liquid  zone  is  increased  by  the  presence  of 
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T  (  °C) 


z  (m) 


Fig.  7  . 


Comparison  of  the  temperature  profiles  in  the 
liquid  zone  between  tap  water  and  degassed 
water.  The  flow  rates  are  identical. 
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gas  in  solution  at  atmospheric  pressure  which  is  released 
under  vacuum  to  form  tiny  bubbles  creating  turbulence  in  the 
liquid.  Unlike  the  atmospheric  pressure  case  or  when  using 
degassed  water,  the  temperature  rises  immediately  at  the 
bottom  of  the  tube. 

A  new  set  of  equations  for  the  liquid  zone  using  water 
from  the  tap  was  developed: 


In (T* ) 

=  Z*  (B,  + 

C,  .Z*) 

(49) 

with  B, 

=  0.3851 

-  6739.3/Pei 

(50) 

and  C , 

=  -0.081 

-  1719/Pei  +  3.12  10 7/Pei  2 

(51  ) 

In  this  case  the  Nusselt  Number  is  given  by: 


Pe  ; 

Nuz  =  - (B,  +  2C ,  .Z*)  (52) 

4 


This  equation  permits  the  evaluation  of  the  local  heat 
transfer  coefficient  in  the  liquid  zone  when  gas  bubbles  are 
released  under  vacuum. 

Equation  44  or  49  will  be  used  in  the  testing  of  the 
model  to  determine  when  the  boiling  temperature  is  reached. 


. 
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b .  Boil i ng  zone . 

In  the  Theory  section  we  developed  the  equations 
controlling  the  pressure  drop,  the  void  fraction  and  the 
heat  transfer  in  the  boiling  zone.  Now  we  have  to  solve 
these  equations  to  see  if  they  lead  to  values  consistent 
with  the  experiment. 

Experimentation  in  the  boiling  zone  was  conducted  with 
water  entering  the  tube  at  a  temperature  corresponding  to 
the  boiling  point  in  the  system  (dependent  on  the  vacuum). 

Both  heat  transfer  and  pressure  drop  equations  contain 
the  mass  flow  rate  of  liquid  phase  and  vapor  phase  which  are 
a  function  of  the  position  in  the  tube.  The  only  positions 
where  we  know  these  values  by  experiment  are  at  the  inlet 
and  at  the  oulet  of  the  evaporator  tube. 

As  suggested  by  Wallis  (1969),  we  can  assume  that  the 
quality  increases  linearly  with  the  distance,  therefore  the 
mass  flow  rate  of  the  vapor  increases  linearly  with  the 
distance: 


W2=Z*.W2exit  (53) 


Differentation  with  respect  to  Z*  gives: 


dW  2 

-  =  W  2  e  x  i  t 

dZ* 


(54) 


.  ■ 
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Now  Equation  37  can  be  solved  numerically  using  a  4th  order 
Runge  Kutta  method  to  obtain  local  values  of  the  pressure 
along  the  tube.  Then  the  temperature  in  the  vapor  phase  can 
be  calculated  from  Equation  31. 

As  the  distance  Z*  is  incremented  at  each  step  of  the 
numerical  integration,  a  local  heat  transfet  coefficient  Uz 
is  evaluated  from  Equation  43. 

The  detailed  computer  program  is  given  in  Appendix  B. 

Figure  8  shows  the  approximate  temperature  profile 
obtained  by  numerical  solution  compared  to  the  experimental 
temperatures  from  Table  2.  The  values  of  the  local  heat 
transfer  coefficient  Uz  obtained  from  all  data  (variation  of 
feed  rate,  variation  of  steam  pressure,  variation  of  vacuum: 
Tables  2,  3  and  4)  appeared  to  be  the  same  for  a  given 
position  Z*.  Table  5  show  the  mean  and  the  standard 
deviation  of  the  estimated  values  of  Uz  while  Figure  9 
represents  the  best  least  squares  fit  of  the  average  heat 
transfer  coefficient  versus  Z*  for  a  regression  of  degree  2. 
(square  of  multiple  correlation  coefficient  =  0.9999). 


Uz  =  615.5  -  1  12  Z*  +  41 .55  Z*2  (55) 

Uz  can  be  also  related  to  a  in  the  following  way: 

Uz  =  615.5  W.itt  2  .  K"  1  for  a  <  0.64 


Uz 


805.44  -  295.4  a  for  0.64  <  a  <  1 


56 
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T(  °C) 


Fig. 8.  Temperature  profile  in  the  boiling  zone: 

The  solid  curve  represents  the  numerical  solution. 
The  vertical  lines  represent  the  experimental 
values  as  per  Table  2. 


Feed  rate  Vapor  rate  T.C.2  T.C.3  T.C.4  T.C.5  T.C.6  T.C.7 

kgs'.  103  kg  s  V 1 03  °C  °C  °C  °C  °C  °C 
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Table  2.  Temperature  profiles  in  the  boiling  zone  for  seven  different  feed 
rates.  The  steam  pressure  and  the  vacuum  remain  the  same. 


47 


00 

CO 

CO 

CD 

Is* 

do 

ID 

ID 

ID 

ID 

i-  e 

ID 

LO 

ID 

ID 

ID 

CD 

CO 

CO 

O 

O 

CO 

<->0 

CD 

CD 

S 

S 

S 

K  ° 

LO 

ID 

ID 

ID 

LD 

ID 


O 


CO 

uy 


o 


U) 


Q.O> 

CO 

> 


CO 

CD 

L-  ^ 

°-o 

E 

co 


O 

CO 

LD 

LD 

S 

00 

CO 

00 

00 

ID 

ID 

LD 

LD 

ID 

00 

CM 

O 

CM 

CO 

0) 

cri 

cri 

CD 

ID 

LD 

LD 

ID 

ID 

ID 

ID 

00 

LD 

c- 

O 

o 

o 

o 

o 

CD 

CD 

CD 

CD 

CD 

CM 

CM 

CO 

CM 

d  o 

CM 

CM 

CM 

CM 

CM 

1- c 

CD 

CD 

CD 

CD 

CD 

CD 

rj- 

CO 

cn 

CD 

T— 

CM 

<3- 

CM 

CO 

CO 

CO 

CO 

CD 

CO 

CD 

CM 

CO 

«— 

CO 

CM 

CD 

O 

<— 

CM 

LD 

(D 

■— 

> — 

*— 

c/) 

© 

k_ 

D 

w 

(/> 

(1) 


E 

03 

<D 

</3 


C 

(1) 

i_ 

03 


<1) 

E 

03 

CO 

<D 

.c 


c 

©« 

o  E 

N  © 

?  © 
—  © 
o  u 

-“■D 

03  © 
2  © 
-i—  C4_ 

.5  5 

CO 

©  TJ 

—  c 

o  ro 

o.E 

e  3 

3  © 

■£  > 
TO 

a3  © 
© 


TO 

I- 


48 


o 

CO 

CO 

CO 

CJCJ 

LD 

cm' 

coo 

CM 

H  ° 

LO 

CD 

CD 

CD 

CD 

ID 

CO 

ID 

CO 

CJCJ 

CD 

CO 

crj 

CO 

ID 

h-  ° 

LD 

CO 

CD 

ID 

LD 

CO 

CO 

CM 

CO 

'Ct 

o 

CD 

ID 

CD 

LD 

LD 

ID 

O 

*— 

cri 

ID 

LD 

r< 

ID 

CD 

h' 

CO 

oy 

CO 

CO 

O 

o 

o 

d 

CD 

CM 

co 

CO 

i- 

CD 

CD 

CM 

CO 

O 

LD 

C 1 3 

CD 

OO 

CM 

00 

CO 

CD 

co 

1- 

CD 

CD 

CD  O 
+->  < — 
CO 

(/3 

D  ^ 

a.-* 

to 

> 


o 

CO 

CO 

*— 

co 

o 

CO 

CO 

CM 

CM 

CM 

<r- 

U) 

E 

3 

3 

O 

CO 

> 


C  0) 

a)  x: 

CD 


—  co 

^  E 

i _  CD 


CD 

CD  TO 
C  k- 
O 

N  T3 
CD 
CDO) 
C 

o  ® 


CD  T3 
JC  C 
CO 


.t  CD 


3 

0) 

CO 

CD 


CO 
03 

O  a 

°"E 

CD  CO 

3  +-* 
CO 
CO 

<d  a> 

I- 

CD 


Eo 

ID 

ID 

ID 

co 

03 

3  ir: 

< — 

» — 

co 

CD 

D 

o  to 

CM 

h- 

03 

CM 

03  CL 

CM 

CM 

co 

CO 

> 


03 

_Q 

CO 


49 


Table  5.  Mean  and  standard  deviation  of  the  heat  transfer 
coefficient  in  the  boiling  zone  from 
experiments  as  per  Tables  2,  3  and  4. 
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Best  least  squares  fit  of  the  heat  transfer 
coefficient  in  the  boiling  zone. 
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VI.  TESTING  THE  MODEL 


When  the  evaporator  is  used  under  vacuum  and  the  feed 
introduced  at  a  temperature  lower  than  the  boiling 
temperature,  the  known  parameters  are: 

-initial  conditions  (T  ,  Cp  ,  kj) 

-feed  rate  W. 

-vacuum  in  the  system  Pv . 

-steam  pressure  Ps  . 

-Equations  44  and  49. 

The  parameters  to  be  determined  are  the  following  ones: 
-Zb*  =  distance  at  which  boiling  starts. 

-W2exit  =  vapor  rate  at  the  outlet  of  the  evaporator. 
-Texit  =  temperature  at  the  outlet  of  the  evaporator. 

A.  Distance  at  which  boiling  starts. 

To  determine  the  temperature  profile  in  the  liquid 
zone,  depending  on  the  presence  or  absence  of  gas  released 
under  vacuum,  we  have  to  solve  Equation  12  using  either 
Equation  44  (no  gas)  or  Equation  49  (with  gas)  to  evaluate 
ln(T*)  as  a  function  of  the  initial  Peclet  Number  and 
position. 

1.  No  bubbles  in  the  liquid  zone. 

Equation  12  gives: 


(57) 
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Introduction  of  Equation  44  gives: 


4  fz*  Pe-t 

In ( T* )  =  -  /  -  (B  +  2C.Z*)dZ* 

Pe,  J  Zc*  4 


Integration  with  respect  to  Z*  gives: 


Z  * 

In ( T* )  =  [B.Z*  +  C.Z*2  ]  (58) 

Z  o  * 


In ( T* )  =  -B . Z  0  *  -C.Z0*2  +  B.Z*  +  C.Z*2 


(59) 


The  computer  program  of  Appendix  C  shows  the  case  when 
no  gas  is  released. 

The  coefficient  A  is  evaluated  from  an  average  value  of 
Z0*  for  the  five  lower  flow  rates  as  given  in  Table  1: 

A  =  -  Z  o  * ( B  +  C.Zo*) 

The  values  for  the  two  higher  flow  rates  were  not  used 
because  such  high  flow  rates  are  not  encountered  when  using 
the  evaporator  under  vacuum;  the  temperature  at  the  exit 
being  lower  than  the  boiling  point,  therefore  no  evaporation 
occurs  in  the  tube. 

-From  the  vacuum  in  the  system  we  know  the  temperature 
at  which  boiling  starts  (the  vapor  follows  the  saturation 
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curve ) . 

-The  boiling  point  rise  due  to  the  concentration  of  the 
feed  has  to  be  added  when  not  using  pure  water. 

-Solution  of  Equation  59  for  that  temperature  gives  us 
the  distance  Zb*  at  which  the  boiling  zone  begins. 

2.  Bubbles  released  in  the  liquid  zone. 

Introduction  of  Equation  49  in  Equation  57  gives: 


fz* 

In  ( T*  )  =  /  (B,  +  2C,  .  Z*  )  dZ* 

J  o 


(60) 


Integration  with  respect  to  Z*  gives: 


In ( T* )  =  B, .Z*  +  C, .Z*2  (61 ) 


Solution  of  Equation  61  for  the  boiling  temperature  as 
defined  earlier  gives  us  the  distance  Zb*  at  which  boiling 
starts. 

B.  Temperature  and  vapor  rate  in  the  boiling  zone. 

In  the  Results  section,  we  used  the  experimental  value 
of  W2exit  to  determine  dW2/dZ*.  Now  dW2/dZ*  is  unknown  and 
we  have  to  use  Equation  43  to  determine  the  change  of  vapor 
flow  rate  with  the  position: 
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dW2  7T.D.C.L.U 


(W,.Cp,  +  W2 .Cp2 )n.c .Pn" 1  dP 


dZ  *  X 


(Psn  -  P")- 


(62) 


X 


dZ* 


with  Uz  given  by  Equation  55. 

-starting  at  Z*=Zb*  we  solve  numerically  Equations  37 
and  62  using  a  double  Runge  Kutta  4th  order  method. 

-When  Z*=1  is  reached,  we  know  the  pressure  at  the  exit, 
therefore  the  temperature,  and  we  know  the  vapor  rate 
W  2  e  x  i  t . 

The  computer  program  to  evaluate  these  parameters  is 
given  in  Appendix  C. 

Table  6  shows  the  comparison  of  experimental  results 
with  the  results  obtained  from  the  model  using  Equation  49 
in  the  liquid  zone  and  city  water  as  medium. 

We  can  see  that  the  agreement  is  very  good,  the  maximum 
difference  for  the  distance  at  which  boiling  starts  being 
3.5%  while  the  maximum  difference  on  the  vapor  rate  at  the 
exit  of  the  evaporator  is  4.43%. 

Table  7  shows  results  obtained  for  degassed  water  using 
Equation  44  in  the  liquid  zone. 

Even  when  using  degassed  water,  the  boiling  starts  a 
little  bit  earlier  than  predicted  numerically,  probably 
because  some  gas  is  still  present  in  the  feed  (the  feed  tank 
being  held  at  atmospheric  pressure,  some  gas  went  back  into 
solution  during  the  process). 
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Table  7.  Comparison  between  experimental  values  and  values  predicted  by  the 
model  for  degassed  water  and  different  liquid  solutions. 
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c.  Testing  with  liquid  food. 

The  model  was  tested  with  sucrose  solutions  of 
different  concentrations  and  tomato  juice.  We  chose  tomato 
juice  because  it  is  well  known  in  the  food  industry  as  being 
a  product  difficult  to  evaporate  due  to  the  high  proportion 
of  pulp  it  contains.  The  results  are  shown  in  Table  7. 

The  sucrose  solutions  were  not  degassed  and  the 
difference  in  the  boiling  position  between  experiment  and 
calculation  decreased  when  the  concentration  of  the  feed 
increased.  This  is  due  to  the  viscosity  effect;  the  more 
viscous  the  product  is,  the  more  difficult  it  is  for  the  gas 
to  be  released. 

For  tomato  juice  we  could  not  use  the  rotameter  because 
of  the  pulp  present  in  the  feed  and  the  color.  The  feed  rate 
was  determined  by  recording  the  decrease  in  weight  of  the 
feed  tank  during  the  process. 

Even  with  the  problems  created  by  the  release  of  gas 
bubbles  in  the  liquid  zone  under  vacuum,  we  can  see  from 
Table  7  that  the  differences  between  the  experimental  values 
and  the  computed  values  for  the  distance  at  which  boiling 
starts  and  for  the  vapor  rate  at  the  exit  of  the  evaporator 
is  always  below  10%.  This  is  very  good  considering  that  the 
error  on  the  feed  rate  reading  can  be  3%  and  the  error  on 
the  concentrate  rate  is  between  2%  and  5%  depending  on  the 
rate  at  which  condensate  enters  the  receiver.  The  error  on 
the  vapor  rate,  observed  experimentally  at  the  exit,  is  then 
between  5%  and  8%. 
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VII.  DISCUSSION  OF  RESULTS 


A.  Liquid  zone. 

In  the  development  of  the  equations  to  determine  the 
heat  transfer  in  the  liquid  zone,  the  Peclet  Number  came  out 
naturally  during  the  dimensionalizat ion  of  the  equations. 

All  the  initial  conditions  are  known:  temperature  of 
the  feed,  feed  rate,  heat  capacity  and  thermal  conductivity 
of  the  feed  at  the  inlet  temperature;  but  some  of  these 
quantities  vary  along  the  tube  and  are  therefore  unknown  for 
a  position  Z*  ^  0.  This  is  why  it  is  advantageous  to  develop 
the  equations  using  the  initial  Peclet  Number  (Pei),  for  it 
is  the  only  Peclet  Number  we  can  calculate  and  the  one  which 
will  reflect  the  initial  conditions. 

In  the  liquid  zone,  we  obtained  two  equations 
correlating  the  Nusselt  Number  to  the  initial  Peclet  Number 
and  the  position  Z*.  These  two  equations  are  of  a  similar 
form  but  differ  in  their  coefficients.  This  is  due  to  the 
fact  that  we  used  the  thermal  conductivity  of  liquid  water 
in  both  cases,  while,  when  gas  bubbles  are  present,  the  real 
value  is  somewhere  between  the  thermal  conductivity  of  the 
liquid  and  the  thermal  conductivity  of  the  gas. 

In  Chapter  13  of  Bird  et  al .  (1960),  the  solution  to 

the  Graetz  equation,  as  proposed  by  Leveque ,  applies  to 
laminar  flow  in  pipes.  It  is  as  follows: 
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Nuin  =  1  .  62  (Re  .Pr  .D/L)  0  •  3  3 
with  Nuin  =  hin  .D/k 


In  our  case  the  Nusselt  Number  is  also  a  function  of  the 
Peclet  Number  (Re.Pr).  The  flow  in  the  liquid  zone  being 
laminar,  this  supports  the  validity  of  our  model  because  the 
equation  we  developed  is  similar  to  the  Leveque  equation. 

In  the  liquid  zone  when  the  system  is  under  atmospheric 
pressure  or  when  using  degassed  water,  we  found  that  the 
temperature  at  the  center  does  not  increase  immediately.  All 
the  flow  rates  investigated  correspond  to  the  laminar  flow 
zone  (Re  <  2100).  As  is  well  known,  the  velocity  profile  is 
parabolic  for  laminar  flow  in  tubes.  The  temperature  profile 
for  radial  heat  flow  from  the  wall  of  a  tube  is  laminar  as 
well;  this  explains  why  the  thermocouple  which  is  situated 
at  15  cm  from  the  beginning  of  the  steam  jacket  at  the 
center  of  the  tube  records  almost  the  same  temperature  as 
the  thermocouple  located  before  the  heating  zone.  However, 
after  a  short  distance  of  heating,  natural  convection  occurs 
and  the  velocity  and  temperature  profiles  are  modified.  In 
Figure  10  the  solid  curve  represents  the  parabolic  velocity 
profile  encountered  in  laminar  flow,  while  the  short  dotted 
curve  represents  the  theoretical  velocity  profile  which 
would  be  encountered  in  the  heated  tube  filled  with  water 
(with  a  zero  flow  rate)  subject  to  natural  convection.  The 
long  dotted  curve  represents  the  resultant  velocity  profile 
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Temperature  profiles  in  the  liquid  zone. 

Solid  curve  =  Laminar  flow  profile. 

Short  dotted  curve  =  natural  convection  profile. 
Long  dotted  curve  =  resultant  profile. 


Fig.  10. 
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which  is  supposed  to  occur  in  the  tube  after  a  certain 
distance  of  heating.  It  can  be  seen  that  this  profile,  which 
is  the  sum  of  the  two  previous  ones,  is  no  longer  parabolic 
and  is  close  to  the  velocity  profile  encountered  in 
turbulent  flow. 

When  working  under  vacuum,  the  release  of  the  dissolved 
gases  prior  to  the  heating  zone  creates  a  turbulent  flow 
velocity  profile  before  the  liquid  enters  the  evaporator 
tube,  and  then  the  temperature  at  the  center  rises 
immediately  as  confirmed  by  Figure  7.  This  phenomenon  is  not 
as  critical  for  solutions  of  high  viscosity,  like  sugar 
solutions  and  tomato  juice,  as  for  water  because  the  rising 
velocity  of  the  bubbles  is  lowered  by  the  viscous  drag  and 
for  highly  concentrated  sucrose  solutions  and  tomato  juice 
the  bubbles  travel  almost  at  the  same  speed  as  the  liquid, 
reducing  their  agitation  effect  on  the  velocity  profile. 

In  the  Results  section,  the  model  was  validated  using 
water.  The  dimensionalizat ion  of  the  equations  was  done  to 
take  into  account  the  influence  of  the  change  in  physical 
properties  which  will  appear  when  using  other  liquid 
products  (viscosity,  heat  capacity,  density  and  thermal 
conductivity).  This  is  why  the  model  gave  good  agreement 
with  the  experimental  values  when  tested  with  sucrose 
solutions  and  tomato  juice. 
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b .  Boil  i ng  zone . 

In  the  boiling  zone  a  first  attempt  at  modeling  was 
made  using  the  equations  proposed  by  Wallis  (1969)  which  are 
concerned  with  separated  flow  with  phase  change.  The  model 
did  not  work  regardless  of  the  repartition  of  the  force  for 
momentum  increase  due  to  phase  change  between  Stream  1  and 
Stream  2.  The  obtained  pressure  drop  was  always  too  high  and 
the  corresponding  temperature  at  the  exit  of  the  evaporator 
too  low.  Even  in  some  cases  the  predicted  exit  temperature 
was  negative  (below  freezing). 

In  the  equations  proposed  by  Wallis,  the  drag  force 
exerted  by  the  tube  wall  on  the  vapor  and  on  the  liquid  and 
the  drag  force  between  the  two  components  (liquid  and  vapor) 
have  to  be  evaluated  using  the  Lockhart  and  Martinelli 
(1949)  correlation  which  was  established  for  isothermal 
two-phase,  two-component  flow  without  phase  change.  In  our 
case  these  two  hypotheses  cannot  be  made:  the  flow  is  not 
isothermal  and  the  phase  change  is  very  important. 

The  equations  proposed  by  Malnes  (1975)  were 
investigated  together  with  the  assumption  made  by  Levy 
(1960)  that,  when  the  rate  of  evaporation  is  rapid  enough, 
friction  and  body  terms  can  be  neglected  as  compared  to 
inertia  and  phase  change. 

We  made  the  same  assumption  as  Malnes  that  the  slip 
ratio  stays  constant. 

The  slip  ratio,  S,  is  given  by:  S=V2/V, 
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The  introduction  of  V,  and  V2  as  per  Equations  20  and 
21  gives : 


p(  G  2  (1  -  a) 

p  2  G  !  CL 


The  mass  fluxes  are  related  to  the  mass  flow  rates  as 
follows : 


G,  =  W ,  /  ( 7T.  D 2 /4  )  and  G2  =  W2/(7T.D2/4) 


We  have  now: 


Pi  W2  (1  -  a) 
p  2  W  i  (X 


We  can  see  that  an  increase  in  W2  is  compensated  by  an 
increase  in  a,  while  the  decrease  in  the  value  of  W,  is 
compensated  by  a  decrease  in  (1  -  a).  Then  the  expression 
W 2 /W , . ( 1 -a ) /a  remains  constant. 

The  slip  ratio  is  then  proportional  to  p,/p2  which  can 
be  assumed  constant  for  a  given  experiment. 
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For  this  model  the  void  fraction  had  to  be  determined 
from  another  equation.  The  Darmouth  correlation  modified  by 
Turner  (1966)  was  appropriate  because  it  is  given  for 
vertical  upward  cocurrent  annular  flow.  It  is  developed 
assuming  that  the  regime  in  the  annular  flow  presents 
viscous  and  surface-tension  forces  which  are  small,  and  the 
fluid  dynamics  are  governed  by  a  balance  between  hydrostatic 
forces  and  inertia  forces  in  the  gas  and  the  liquid.  This  is 
compatible  with  the  assumption  made  by  Levy  (1960)  that 
inertia  and  phase  change  terms  are  usually  much  larger  than 
friction  and  body  terms  when  rapid  evaporation  occurs  in  the 
tube . 

This  choice  of  calculation  for  a  was  further  supported 
by  the  fact  that  all  our  experimental  values  for  the 
dimensionless  volumetric  fluxes  j,*  and  j2*  were  within  the 
graph  given  by  Wallis  (1969)  representing  the  range  for 
which  the  equation  proposed  by  Turner  (1966)  is  valid. 

c.  About  the  computer  programs. 

In  the  numerical  evaluation  of  the  pressure  in  the 
boiling  zone  (see  Appendix  B)  the  step  size  used  in  the 
Runge-Kutta  method  was  set  at  1/20.  This  choice  is  not  an 
arbitrary  one.  A  series  of  tests  were  made  on  the  same  set 
of  data  for  different  step  sizes  of  1/5,  1/10,  1/20  and 
1/50.  The  numerical  values  of  pressure,  heat  transfer 
coefficient  and  void  fraction  at  the  exit  of  the  evaporator 
( Z *=  1  )  were  compared.  The  difference  between  values  obtained 
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for  a  step  size  of  1/5  and  values  obtained  for  a  step  size 
of  1/10  was  about  10%;  the  difference  between  a  step  size  of 
1/10  and  a  step  size  of  1/20  was  2%;  while  the  difference 
between  a  step  size  of  1/20  and  a  step  size  of  1/50  was  less 
than  0.3%.  A  difference  of  0.3%  being  a  lot  smaller  than  the 
precision  we  have  on  experimental  data,  a  step  size  of  1/20 
was  chosen.  A  step  size  of  1/50  could  have  been  used,  but  it 
would  have  required  more  computing  time  and  greater  expense 
for  a  gain  in  precision  which  is  not  significant.  As  stated 
earlier  (Testing  the  model  section),  the  error  on  the  vapor 
rate  observed  exper imen taly  at  the  exit  of  the  evaporator  is 
between  5%  and  8%. 

In  Appendix  C  a  step  size  of  (l-Zb*)/20  is  used.  This 
is  smaller  than  1/20  and  we  are  sure  that  the  precision  of 
the  numerical  solution  is  at  least  as  good  as  the  one 
obtained  in  Appendix  B. 
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VIII.  APPLICATION  TO  AN  INDUSTRIAL  SIZE  EVAPORATOR 

Data  for  an  industrial  size  evaporator  were  provided  by  Dr. 
Merson  of  the  Department  of  Food  Science  and  Technology, 
University  of  California,  Davis,  California.  They  are 
summarized  in  Figure  11. 

a.  Description  of  the  evaporator. 

A  schematic  diagram  of  the  evaporator  to  be  studied  is 
given  in  Figure  11. 

This  evaporator  is  used  to  concentrate  pineapple  juice 
from  13.5°Brix  to  62°Brix  at  a  rate  of  approximately 
5,000  kg/hr . 

As  can  be  seen  on  the  diagram,  the  evaporator  is 
composed  of  3  consecutive  sections,  each  section  having  more 
tubes  than  the  preceeding  one  to  permit  expansion  of  the 
vapor  and  ensure  a  climbing  film  pattern  during  the  entire 
process . 

The  only  known  parameters  are  the  inlet  conditions,  the 
cutlet  conditions  and  the  characteristics  of  the  evaporator 
(tube  diameter,  length,  thickness  and  number). 

The  feed  is  preheated  before  entering  the 
evaporator  (115°C). 

The  pressure  at  the  outlet  is  known  (495mm  Hg ) .  The 
pressure  at  the  inlet  of  the  evaporator  is  unknown,  and  an 
appropriate  value  can  be  estimated  by  analogy  with  our 
single  tube  glass  evaporator.  The  average  pressure  drop  in 
our  glass  evaporator  is  47.6mm  Hg  for  the  2.67m  length.  The 
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Fig.  1 1 . 


Schematic  diagram  of  the  industrial  evaporator. 
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added  tube  length  of  the  industrial  evaporator  is  6.7m;  the 
expected  pressure  drop  will  be  close  to  119.5mm  Hg ,  which 
added  to  the  outlet  pressure,  gives  an  approximate  inlet 
pressure  of  614.5mm  Hg .  This  corresponds  to  an  inlet 
temperature  of  94°C.  Because  of  that,  flashing  occurs  in  the 
conical  section  preceding  the  first  set  of  tubes,  so  a 
mixture  of  vapor  and  liquid  reaches  the  first  section  of  the 
evaporator  and  two  phase  flow  is  present  all  along  the  tubes 
in  each  section. 

b.  Method  of  cal cul at  ion . 

In  order  to  apply  our  model  to  that  evaporator  we  have 
to  consider  each  section  separately. 

We  assume  that  the  flow  is  equally  distributed  among 
the  tubes  in  each  section.  We  divide  the  total  flow  by  the 
number  of  tubes  present  in  that  section. 

We  apply  the  model  for  one  tube  (after  modification  for 
the  character i st ics  of  that  tube:  length,  diameter  and  heat 
transfer  coefficient).  This  permits  one  to  evaluate  the 
pressure  at  the  top  of  that  section,  while  the  vapor  and 
concentrate  rates  are  determined  by  multiplying  the  results 
for  one  tube  by  the  number  of  tubes  in  that  given  section. 

The  composition  at  the  exit  of  that  section  being  known 
(void  fraction  and  density  of  both  phases),  the  pressure 
drop  in  the  space  between  tube  sheets  from  one  section  to 
the  next  one  can  be  evaluated.  Assuming  homogenous  two-phase 
flow,  the  pressure  drop  between  sections  is  given  by: 


■ 
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AP  =  H  [  a  .  p  2  +  (1  "  a)p,  3  .g 


(63) 


with  H  =  0.153m 


1.  Evaluation  of  the  amount  of  flashing. 

The  amount  of  flashing  prior  to  the  first  section  of 
the  evaporator  is  not  known,  and  the  composition  of  the  feed 
(ratio  of  vapor  to  liquid)  and  the  temperature  at  the 
beginning  of  the  evaporation  process  are  unknown. 

The  only  parameter  available  is  the  pressure  at  the 
exit  of  the  evaporator,  therefore  the  temperature  is  known 
(the  vapor  phase  is  assumed  to  follow  the  saturation  curve). 

An  iterative  evaluation  converging  on  the  outlet 
temperature  has  to  be  used.  It  works  as  follow: 


V  W 


The  initial  conditions  prior  to  the 
inlet  of  evaporator  are  known: 

Temperature,  feed  rate,  degree  Brix. 
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Cone  section 


Section  1 


Equation  63 — >AP 


FLASHING  (unknown) 

I  * . 

W,  ,  W 2  ,  P  <  —  Guess  on  temperature  <- 


Equations  37  and  62 


W ,  ,  W  2  ,  P 


v 

Equations  37  and  62 

Y 

W, ,  w2 ,  P 


Y 

Equations  37  and  62 


Y 

W, ,  W2 ,  P 

- - - -  i 

°Brix  T  estimated  (Test.) 

Convergence  test: 

-  if  Test.  ^  Texp. 


-  if  Test.  =  Texp. 


RESULTS 


W,  ,  W2 ,  0 Br i x 


Exit  temperature 
i s  known  ( Texp . ) 


Sect  ion  2 


Equation  63 — >AP 


Section  3 


EXIT 
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2.  Evaluation  of  the  heat  transfer  coefficient. 

When  the  model  was  set  up  for  the  single  tube  glass 
evaporator,  an  overall  heat  transfer  coefficient  between  the 
steam  outside  the  tube  and  the  vapor  inside  the  tube  was 
evaluated.  The  resistance  to  heat  transfer  is  composed  of 
three  resistances  in  series:  one  through  the  film  of 
condensed  steam  on  the  steam  side,  one  through  the  wall  of 
the  tube  itself  and  one  through  the  film  of  liquid  climbing 
the  inside  of  the  evaporator  tube. 

We  have  to  use  the  values  obtained  for  the  single  tube 
glass  evaporator  to  determine  the  heat  transfer  coefficient 
of  the  industrial  evaporator. 

If  we  know  the  thickness  of  the  wall  and  the  material 
it  is  made  of,  the  resistance  through  the  wall  can  be 
evaluated  but  the  resistances  on  the  steam  side  and  on  the 
film  side  are  not  directly  available. 

The  diameter  of  the  tubes  being  much  bigger  than  the 
wall  thickness,  we  evaluate  the  heat  transfer  coefficient 
using  the  flat  wall  equation,  which  is  as  follows: 

1/U  =  1/hs  +  Ax/k  +  1/hf  (64) 


For  the  glass  single  tube  evaporator  we  have: 


1/UA 


1/hs  +  A  x A /k A  +  1/hf 


(65) 
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with  axa  =  1.53  1  0 "  3m  and  kA  =  1.163  W.nr’.K'1 


For  the  stainless  steel  evaporator  we  have: 


1/UB  =  1/hs  +  A  x  B /kB  +  1/hf 


(66) 


with  AxB =  1.27  10'3m  and  kB  =  15  W.nr’.K"1 


Let  us  assume  that  the  heat  transfer  coefficients  on  the 
steam  side  are  similar  in  both  case  and  that  the  heat 
transfer  through  the  film  is  the  same  for  a  given  void 
fraction  a.  Elimination  of  those  two  factors  between 
Equations  65  and  66  gives: 


1/UB  =  1/UA  +  AxB/kB  -  AxA/kA  (67) 

UA  is  known  for  different  values  of  the  void  fraction  a. 
(refer  to  Table  5).  Then  1/Ub  and  therefore  Ub  can  be 
evaluated  as  a  function  of  a.  UB  is  given  by  the  following 
equa t i on : 


[i  t ; 
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UB  =  2539  W.m~  2 . K" 1  for  a  <  0.619 

UB  =  4657.1  -  3423.7  a  for  0.619  <  a  <  1 


(68) 


The  square  of  the  multiple  correlation  coefficients  for  the 
first  degree  best  least  squares  fit  is  0.947. 

3.  Results. 

The  program  converged  after  5  iterations.  The  values  of 
the  concentrate  rate,  temperature  and  °Brix  at  the  exit  of 
each  section  are  given  in  Table  8. 

The  values  for  Section  3  correspond  to  a  °Brix  of 
62.24,  which  is  only  0.4%  higher  than  the  experimental 
62  0  Br i x . 

Looking  at  the  vapor  rate  difference  between  the  model  value 
and  the  experimental  value,  we  have  an  experimental  rate  of 
3917  kg.hr'1  and  a  calculated  rate  of  3919  kg.hr'1.  The 
difference  is  only  0.05%. 
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CONCLUSIONS 

The  results  of  this  investigation  indicate  that  it  is 
possible  to  predict  and  control  the  performance  of  a 
climbing  film  evaporator,  assuming  the  characteristics  of 
the  material  to  be  concentrated  are  known. 

The  heat  transfer  in  the  lower  liquid  zone  is  the 
easiest  to  solve  because  it  is  controlled  by  very  well  known 
parameters  (in  this  zone  the  evaporator  works  like  a  heat 
exchanger).  It  was  more  difficult  to  solve  the  problem  for 
the  two-phase  flow  region  where  no  successful  theory  was 
previously  developed  for  this  kind  of  evaporator. 

In  the  boiling  zone  where  flashing  occurs,  our  theory 
is  developed  from  the  equations  given  by  Malnes  (1975),  who 
worked  on  superheated  liquids  encountered  in  water  cooled 
nuclear  reactors.  His  study  was  done  for  very  high  pressure 
conditions  in  the  system,  while  the  present  work  is  done 
under  vacuum  conditions,  making  these  equations  universal 
for  evaporation  in  vertical  tubes  when  flashing  is  present. 

This  is  further  confirmed  by  the  ability  of  the  model 
to  predict  the  performance  of  a  big  industrial  evaporator. 

The  model  can  be  used  to  optimize  the  performance  of  a 
given  evaporator  by  determination  of  the  best  parameters 
(feed  temperature,  feed  rate,  vacuum  and  steam  pressure)  for 
the  concentration  of  a  product  from  a  given  initial 
concentration  to  a  desired  final  concentration. 


esw  ^  *******  to  *aarig  ■  *d$ 
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The  model  can  also  be  very  useful  for  the  design  of  an 
evaporator  to  accomplish  a  specific  evaporation  during  a 
given  process.  The  thermal  properties  and  the  amount  of 
product  to  be  concentrated  being  given,  the  model  can  be 
adapted  to  determine  the  size  of  the  evaporator  (number  of 
tubes,  number  of  sections  and  tube  length)  to  be 
manufactured . 
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APPENDIX  A. 
Data . 


A.  Temperature  vs  pressure  for  the  vapor  phase. 

The  values  of  temperature  vs  pressure  as  given  by 
Loncin  and  Merson  (1979)  were  used  to  derive  equation  31. 
The  temperatures  were  converted  to  absolute  temperatures 
(K)  . 

The  first  degree  best  least  square  fit  of  ln(T)  vs 
ln(P)  is  as  follow: 

In  ( T )  =  5.11  659  +  0.06978  ln(p) 

The  square  of  multiple  correlation  coefficients  is:  0.9992. 
After  exponentiation  we  obtain: 


T  =  166.766  P°’C6,,B 

8.  Latent  heat  of  vapor izat ion . 

The  values  of  enthalpies  of  vaporisation  as  given  by 
Loncin  and  Merson  (1979)  were  used  after  conversion  to 


. 


. 
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absolute  temperature. 

A  first  degree  best  least  square  fit  gives: 

X  =  3135665.7  -  2345.714  T 

The  square  of  multiple  correlation  coefficients  is  0.9995. 

C .  Boiling  po i nt  r i se  of  sugar  so  1 ut ions . 

The  molecular  elevation  of  the  boiling  point  as  given 
by  Weast  (1974)  is  used.  For  water  it  is  equal  to  0.512°C 
per  g-mol/1.  The  molecular  weight  of  sucrose  is  342.3. 


wt/%wt 

g-mol/1 

T  (  °C) 

17 

0.530 

0.27 

23 

0.736 

0.38 

33 

1.101 

0.56 

40 

1  .375 

0.70 

D .  T hernia  1  proper  t  i  es . 

The  tubes  of  the  industrial  evaporator  are  made  of  18 
gauge  stainless  steel. 

Loncin  and  Merson  (1979)  give  a  thermal  conductivity 
for  stainless  steel  of  15  W.irr'.K'1. 


■ 
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The  thickness  for  18  gauge  U.S.  standard  plate  as  given 
by  Weast  (1974)  is  1.27  10~3m. 

Weast  gives  a  value  of  thermal  conductivity  for 
borosilicate  pyrex  type  glass  of  1.163  W.irr’.K"1. 

The  thickness  of  the  tube  of  our  evaporator  was  not 
directly  available.  The  outside  diameter  was  accessible  for 
measurement.  The  inside  diameter  was  evaluated  by  measuring 
the  volume  of  a  given  height  of  water:  this  volume  was 
divided  by  the  height  to  obtain  the  inside  cross  section 
which  was  used  to  evaluate  the  inside  diameter.  The 
thickness  is  the  difference  between  the  outside  radius  and 
the  inside  radius.  It  was  found  to  be  1.53  10'3m. 

Loncin  and  Merson  (1979)  give  a  value  of 
2,030  J-kg^.K*1  for  the  heat  capacity  of  the  vapor  at 
100°C.  This  is  half  the  value  of  the  specific  heat  of  liquid 
water  at  the  same  temperature.  When  the  concentration  of  the 
liquid  phase  increases,  its  specific  heat  decreases.  Because 
the  heat  absorbed  by  the  liquid  and  the  vapor  phases  is 
small  compared  to  the  heat  of  vaporisation,  in  the  computer 
programs  we  consider  that  the  specific  heat  of  the  liquid  is 
identical  to  that  of  water  while  the  specific  heat  of  the 
vapor  is  taken  as  half  of  this  value. 
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APPENDIX  B. 

Computer  program  for  the  boiling  zone. 


This  computer  program  evaluates  the  pressure  drop  and  the 
local  heat  transfer  coefficient  in  the  boiling  zone  when 
boiling  occurs  all  along  the  tube. 

This  program  is  composed  of  a  main  program  and  two 
subrout ines . 


, 
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C  PRESSURE  DROP  IN  THE  CLIMBING  FILM  EVAPORATOR 

C  IN  THE  BOILING  ZONE. 

COMMON  W 1 , W2 , R 1 , R2 , R , DM2 , T 
C 

C  TO  EVALUATE  THE  PRESSURE  DROP,  THE  VOID  FRACTION 

C  AND  THE  LOCAL  HEAT  TRANSFER  COEFFICIENT  WHEN 

C  FILM  BOILING  IS  PRESENT  ALL  ALONG  THE  TUBE. 

C 

REAL  K 1 ,K2 ,K3 ,K4 
C 

C  K1  TO  R4  ARE  THE  RUNGE  KUTTA  COEFFICIENTS  TO 

C  EVALUATE  DP/DZ* 

C 

READ  (  1  ,  3  0  0  ) K 
300  FORMAT  (12) 

NK=  0 

3  IF(NK.EQ.K)  GO  TO  1 

NK=NK+ 1 

READ  ( 1 , 1 00) PS , DM, DM2 
C 

C  PS=STEAM  PRESSURE  (PSI),  DM=FEED  RATE  (GR/MIN), 

C  DM2 = VAPOR  RATE  AT  THE  EXIT  OF  THE  EVAPORATOR. 

C 

100  FORMAT  (3E12.5) 

READ  (1,110) DPV , X 
C 

c  DPV=VACUUM  READING  (CM  OF  MERCURY), 


. 
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C  X=DI STANCE  AT  WHICH  BOILING  STARTS:  X=0. 

C 

110  FORMAT  (2E12.5) 

R 1  =  1  000 . 

R=0 . 0  1  5 
DM=DM/ 60000 
DM2 = DM2/ 60000 
WRITE  (6, 1 20) PS ,DM,DPV 
120  FORMAT  (/,5X,'ST.  PRESS . = ' , F5 . 2 , 5X , 

* ' FLOW  RATE='  ,E 12.5, 5X ,  ’VAC.  READI NG= '  ,  F7 . 2 ) 

WRITE  (6,130) 

130  FORMAT  ( / , 5X ,  ' Z  = '  ,  1 5X ,  ' TEMP= '  ,  1 2X ,  ' H . T .  COEF.=’,6X, 
* ’ ALPHA= ’ ,/) 

PS= ( 14 . 7+PS) *6890 
N=  2  0 

Y= ( 76 .-DPV)/1 00 . *13600. *9. 81 
Z  =  X 

ALPHA  =  0 . 

TEMP= 166. 766 *Y**0. 06978-273. 

H=( 1 .-Z)/N 
2  W2=DM2  *  Z 

R2  =  Y*  *  0 . 93  0  22/7  7  04  5.8  92 
W1=DM-W2 
P= Y 
C 

C  THE  VALUE  OF  DP/DZ*  IS  EVALUATE  IN 

c  SUBROUTINE  DPM . 


m  '  '  -  i  «fl 
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C  THE  LOCAL  HEAT  TRANSFER  COEFFICIENT  IS 

C  IN  SUBROUTINE  COEF  AFTER  COMPLETION  OF 

C  KUTTA  EVALUATION  OF  DP/DZ*. 

C 

CALL  DPM ( P , DP , U ) 

CALL  COEF (Z, PS, P, TEMP, DP, UZ ) 

WRITE  (6 , 200 ) Z ,TEMP,UZ , ALPHA 
I F ( Z . GT .0.99)  GO  TO  3 
K 1 =H*DP 

R2=(Y+K1/2. )** 0.93022/7 7 04 5. 8 92 
T=  Z+H/2 . 

W2=DM2  *T 
W1 =DM-W2 
P=Y+K1/2 . 

CALL  DPM ( P , DP , U ) 

K2=H*DP 

R2  = ( Y  +  K2/2 . )**0. 93  0  22/7  7  04  5.8  92 
P=Y+K2/2 . 

CALL  DPM ( P , DP , U ) 

K3=H*DP 

R2= ( Y+K3 ) **0.93022/77045.892 

T=Z+H 

W2  =DM2  *T 

W1=DM-W2 

P=Y+K3 

CALL  DPM ( P , DP , U ) 


EVALUATED 
THE  RUNGE 


K4=H*DP 
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Y=Y+(K1+2.*K2+2.*K3+K4)/6. 

Z  =  T 

ALPHA=U 

TEMP= 1  66 . 7 66*Y**0. 06978-273. 
200  FORMAT  ( 4 ( 5X , E 1 2 . 5  ) ) 

GO  TO  2 

1  CONTINUE 

STOP 
END 
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SUBROUTINE  DPM(P,DP,U) 

C 

C  TO  DETERMINE  THE  PRESSURE  DROP  DP/DZ*  IN  THE 

C  BOILING  ZONE. 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , DM2 , T 
REAL  J1 , J2 

IF  (W2.EQ.0. )  GO  TO  1 
C 

C  EVALUATION  OF  THE  VOID  FRACTION  U  BY  SOLVING 

C  EQUATION  38. 

C 

A=  3 .  14  1 6*R*  *  2 
C= 1 . 20736E-05 
A  1  =  3 .  1**2 

D=  9 . 8 1 *0 . 03* (R1-R2 ) 

D=SQRT (D ) 

J 1 =W1 /SQRT ( R 1 )/A/D 
J2=W2/SQRT(R2 )/A/D 
B 1  =  3 . 1* ( J2+J1-1  . ) 

C 1 =- J  1 

DELT=B 1 **2-4 . *A1*C 1 

X= ( ~B 1 +SQRT ( DELT ) )/2./A1 

u=i . -x 

S=W2*R1/W1/R2* ( 1 . -U)/U 
V1=W1/R1/A/( 1 ,-U) 

V2=W2/R2/A/U 


. 
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RO=U*R2*S**2+ ( 1 .-U) *R1 
Z  =  R2  *U  +  R 1  * ( 1 . -U) 

D=V 1  *  *  2  *R0*U*C/R2/P*  *  0 . 06  978 

D=  1  . -D 

DMV=DM2 

B= ( RO* (R1~S*R2)/S/R1 /R2+S- 1 . ) *DMV*V1/A 
B=B  +  9 . 81*2. 67  *  Z 
DP=-B/D 
GO  TO  2 

1  DP=-9 . 8 1 *2 . 67*R1 

2  CONTINUE 
RETURN 
END 
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SUBROUTINE  COEF ( Z , PS , P , TEMP , DP , UZ ) 

C 

C  TO  DETERMINE  THE  LOCAL  HEAT  TRANSFER  COEFFICIENT 

C  IN  THE  BOILING  ZONE  (UZ). 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , DM2 , T 
REAL  LAM 
CP2=4 1 80 
CP  1 =2000 
X  =  TEMP  +  2  7  3  . 

LAM=3 135665. 7-2345. 714*X 
DMV=DM2 

A= (W2*CP2+W1 *CP1 )*0 . 06978*DP/P**0 . 93022 
A=A+DMV*LAM/166. 766 

B=3  .  14  16*2. *R* 2  . 67* ( PS*  *  0 . 06978~P**0. 06978) 

UZ  =a/b 
RETURN 


END 
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APPENDIX  C. 

Computer  program  for  the  model. 


This  program  evaluates  the  performance  of  the  evaporator  in 
mormal  use  given  the  inlet  conditions. 

The  main  program  evaluates  the  distance  at  which 
boiling  starts  and  calls  the  subroutine  VAP  which  evaluates 
the  pressure  and  the  vapor  rate  in  the  boiling  zone  through 
two  other  subroutines:  DPM  and  DVAP. 


. 


•  *  - 
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C 

C 

C 

C 

c 

C 

c 

C 

c 

C 

c 

c 


3 

2  1  0 


220 

1  1  0 


230 

130 


CLIMBING  FILM  EVAPORATOR  PROJECT 
OPERATION  OF  THE  EVAPORATOR  IN  NORMAL  USE. 
INPUT:  STEAM  PRESSURE,  VACUUM  READING. 

INPUT:  FEED  TEMPERATURE,  FEED  RATE 
IN  GRAMS  PER  MIN. 

INPUT:  THERMAL  DIFFUSIVITY,  DENSITY, 

BOILING  POINT  RAISE. 

OUTPUT:  DIMENSIONLESS  LENGTH  WHERE 
BOILING  OCCURS. 

OUTPUT:  VAPOR  RATE  AT  EXIT,  VOID  FRACTION, 
TEMPERATURE 

COMMON  W1 ,W2,R1 , R2 , R , ZB , DW2 , Y , Z 
R= 0 . 0  1  5 
WRITE  (6,210) 

FORMAT  (/,5X,’ STEAM  PRESSURE=?  VACUUM  READI NG=  ? '  ) 

READ  ( 5 , 1 1 0 ) PS , PV 
WRITE  (6,220) 

FORMAT  ( 5X , ' FEED  TEMPERATURE* ?  FEED  RATE=? ' ) 

READ  (5,110) TI , DM 
FORMAT ( 2E1 2 .5) 

DM=DM/ 60000 
WRITE  (6,230) 

FORMAT  ( 5X ,  ’ T . D . =  ?  DENSITY*?  B . P . RAI SE=  ? '  ) 

READ  ( 5 , 1 30 ) TD , R 1 , BPR 
FORMAT  (3E12.5) 

WRITE  (6,240) 


\ 
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240  FORMAT  (//) 

C 

C  CALCULATION  OF  THE  DISTANCE  AT  WHICH  BOILING 

C  STARTS  (ZB)  FROM  EQUATIONS  ESTABLISHED  FOR  THE 

C  LIQUID  ZONE. 

C 

REPR=4 . *DM/R 1 /3 . 1416/0. 03/TD 
B=-0 . 4583+ ( 9444 . 7+5 . 2 1 7E+0 6/REPR ) /REPR 
C=0.2384+(754.6-1. 345E+07/REPR) /REPR 
Z 0=  1  3 . /2 67  . 

A=-b*Z0-C*Z0**2 

PV= (76.-PV)/l00.*13600.*9.81 

TV= 1 66 . 766*PV**0 . 06978-273. 

TV=TV+BPR 

PS=( 1 4 . 7  +  PS ) *6890 . 

TS= 1 66 . 766*PS**0 . 06978-273 . 

TAU= (TS-TI )/(TS-TV) 

TAU=ALOG ( TAU ) 

DELT=B**2+4 .  *C* (TAU- A) 

ZB= ( -B+SQRT ( DELT ) )/2 ./C 
BOI L=  ZB*  2  67  . 

WRITE  (  6 , 2 0 0  ) TI , DM 
WRITE  (6,250) BOIL 

200  FORMAT  (/,5X,'FEED  TEMP= ’ , F6 . 2 , 5X , ' FEED  RATE= ’ , E 1 2 . 5 ) 
250  FORMAT  ( / ,  1  OX ,  ' BOI LI NG  STARTS  AT  Z=\F7.2) 

IF  ( ZB.GE. 1 . )  GO  TO  1 

C 


. 
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C  THE  PRESSURE  AND  THE  VAPOR  RATE  IN  THE  BOILING 

C  ZONE  ARE  CALCULATED  IN  SUBROUTINE  VAP  FROM 

C  Z*  =  ZB*  TO  Z*=  1  . 

C 

CALL  VAP ( DM , PS , PV , TV , U ) 

TEMP= 1 66. 766*Y**0. 06978-273. 

WRITE  ( 6 , 300 ) W2 , U , TEMP 
300  FORMAT  ( / ,  1  OX ,  ' VAPOR  RATE= '  , E 1 2 . 5 , 5X , 

* ' ALPHA= ' ,  E  1  2 . 5 , 5X  , ' TEMP= ' , F5 . 2 , /) 

GO  TO  2 

1  WRITE  (6,400) 

400  FORMAT  (/,3X,'Z3  BIGGER  THAN  ONE,  TRY  NEW  GUESS.’,/) 

2  WRITE  (6,450) 

450  FORMAT  (5X,’SIGN=?  IF  SIGN=0.STOP,  SI GN= 1 . CONTINUE /) 
READ  ( 5  ,  1 20 ) SIGN 
120  FORMAT  (El  2. 5) 

IF  (SIGN.EQ. 1 . )  GO  TO  3 
STOP 


END 


. 
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SUBROUTINE  VAP ( DM , PS , PV , TV , U ) 

C 

C  TO  DETERMINE  THE  PRESSURE  AND  THE  VAPOR  RATE  IN 

C  THE  BOILING  ZONE  USING  A  DOUBLE  RUNGE  KUTTA  4TH 

C  ORDER  METHOD. 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , ZB , DW2 , Y , Z 
REAL  K1 ,K2 ,K3 ,K4 ,L1 ,L2 ,L3 ,L4 
C 

C  K1  TO  K4  ARE  THE  RUNGE  KUTTA  COEFFICIENTS  TO 

C  SOLVE  EQUATION  37. 

C  LI  TO  L4  ARE  THE  RUNGE  KUTTA  COEFFICIENTS  TO 

C  SOLVE  EQUATION  62 

C 

N=  2  0 

H= ( 1 . -ZB)/N 
Z=  ZB 

wv=o. 

Y  =  PV 

3  W2=WV 

P=Y 

W1=DM-W2 

R2  =  P*  *  0 . 93022/77045.892 

C 

C  DP/DZ *  IS  EVALUATED  IN  SUBROUTINE  DPM 

C  DW2 /DZ *  IS  EVALUATED  IN  SUBROUTINE  DVAP 

C 


- 


CALL  DPM ( P , DP , U ) 

CALL  DVAP ( PS , P , DP ) 

K  1  =H*DP 
L 1 =H*DW2 
Z=Z+H/2 . 

W2=WV+L 1/2. 

P=Y+K1/2 . 

R2=P**0. 93022/77045.892 

W1=DM-W2 

CALL  DPM ( P , DP , U ) 

CALL  DVAP ( PS , P , DP ) 
K2=H*DP 
L2=H*DW2 
P=Y+K2/2 . 

W2=WV+L2 / 2 . 

R2=P**0 . 93022/77045.892 

W1=DM-W2 

CALL  DPM ( P , DP , U ) 

CALL  DVAP ( PS , P , DP ) 
K3=H*DP 
L3=H*DW2 
Z  =  Z+H/2 . 

P=Y+K3 

W2=WV+L3 

R2=P**0. 93022/77045.892 

W1=DM-W2 

CALL  DPM ( P , DP , U ) 


■ 

. 
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CALL  DVAP ( PS , P , DP ) 

K4=H*DP 

L4=H*DW2 

Y=Y+ ( K 1 +2 . *K2+2. *K3+K4)/6. 
WV=WV+(L1+2.*L2+2.*L3+L4)/6. 
IF  (Z.GE. 0.999)  GO  TO  2 
GO  TO  3 
2  W2=WV 

RETURN 
END 


SUBROUTINE  DPM(P,DP,U) 


TO  EVALUATE  DP/DZ* 

COMMON  W 1 , W2 , R 1 , R2 , R , ZB , DW2 ,  Y ,  Z 
REAL  J1  ,  J2 

IF  (W2 .EQ. 0 . )  GO  TO  1 
A=3 . 1 4 1 6*R*  *  2 
C=  1  .20736E-05 
A 1  =  3 . 1**2 

D=9.81*0.03*(R1-R2) 

D=SQRT (D ) 

DETERMINATION  OF  THE  VOID  FRACTION  U  BY  SOLVING 
EQUATION  38. 

J 1 =W 1 / SQRT ( R 1 ) /A/D 
J2=W2/SQRT(R2 )/A/D 
B 1  =  3 . 1  * ( J2  +  J1- 1  .  ) 

C 1 =- J 1 

DELT=B 1 **2-4 . *A1 *C1 

W= ( -B 1 +  SQRT ( DELT ) )/2./A1 

U=  1  .  -W 

S=W2*R1/W1/R2* ( 1 .-U)/U 
V1=W1/R1/A/( 1 . -U) 

V2=W2/R2/A/U 
rq=U*R2*S*  *  2+ ( 1 .-U) *R1 


i 
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X=R2*U+R1 * ( 1 ,-U) 

D=V 1  *  *  2  *RO*U*C/R2 /P*  *  0 . 06  978 
D=1  .-D 

B= (RO* ( R 1 -S*R2 )/S/R1/R2+S- 1 . ) *DW2*V1/A 
B=B+9 .81*2. 67*X 
DP=-B/D 
GO  TO  2 

1  DP=-9 . 8 1 *2 . 67*R1 

2  CONTINUE 
RETURN 
END 


103 


SUBROUTINE  DVAP ( PS , P , DP ) 

C 

C  TO  DETERMINE  DW2 /DZ * 

C 

COMMON  W1 ,W2,R1 , R2 , R , ZB , DW2 , Y , Z 
REAL  LAM 

T=  1  66 . 766*P**0 . 06978 
LAM=3 135665. 7-2345. 714*T 
COEF=615. 5-Z*( 112. -41. 5 5* Z) 

C 

C  COEF  =  LOCAL  HEAT  TRANSFER  COEFFICIENT  IN  THE 

C  BOILING  ZONE  AS  DEFINED  IN  EQUATION  55. 

C 

CP=4 180. 

A  =  3  .  1416*0.03*166.766*2.67 
DW2=A*COEF* (PS**0 . 06978~P**0 . 06978 ) 

DW2=DW2- (W2/2 . +W1 ) *CP* 166.766*0. 06978 /P**0 . 93022*DP 

DW2=DW2 /LAM 

RETURN 


END 


’•  3 
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APPENDIX  D. 

Computer  program  for  the  industrial  evaporator. 


This  program  evaluates  the  performance  of  the  industrial 
evaporator  capable  of  concentrating  5000  kg/hr  of  pineapple 
juice  from  13.5°Brix  to  62°Brix. 

This  program  is  composed  of  a  main  program  and  two 
subrout i nes . 
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C  COMPUTER  PROGRAM  TO  EVALUATE  THE  PERFORMANCE 

C  OF  THE  INDUSTRIAL  EVAPORATOR. 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , DW2 ,Y,Z,L,U,V1,V2 
REAL  K 1 ,K2,K3fK4rL1 , L2 , L3 , L4 
H= 1 . /2  0 . 

R=0 . 0095 

PS=(  14.7  +  35. ) *6890 . 

C 

C  INITIAL  GUESS  ON  THE  TEMPERATURE:  88.5  C. 

C  DT  IS  THE  TEST  FOR  CONVERGENCE. 

C  NK  IS  TO  SET  THE  EVAPORATOR  SECTION. 

C  DMF  IS  THE  AMOUNT  OF  INITIAL  FLASHING. 

C 

T=88. 5  +  273. 

DT=0  . 

9  NK=  0 

T=T+DT 

P=(T/166.766)** 14 . 331 

DMF= 1 .39*4000. *(1 1 5 . 5+2 7 3 . -T ) /2 . 2 58E+0 6 
Y  =  P 

W2=DMF/ 66 . 

W1  =  ( 1 . 3  9 -DMF ) /66 . 

C 

C  EVALUATION  OF  THE  INLET  CONDITIONS  FOR 

C  ONE  TUBE  IN  A  GIVEN  SECTION. 

C  R1  IS  THE  DENSITY  OF  THE  LIQUID  PHASE. 


' 

. 
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C  R2  IS  THE  DENSITY  OF  THE  VAPOR  PHASE. 

C 

L= 2  .  1  3 

8  NK=NK+ 1 

IF  (NK.EQ. 1 )GO  TO  1 
IF  (NK.EQ. 2) GO  TO  2 
L=2 . 44 

W1=W1 *  111/156. 

W2=W2* 1 1 1 ./I  56 . 

GO  TO  1 

2  W1=W1*66./1 1 1 . 

W2=W2* 66 . / 1 1 1 . 

1  DM=W 1 +W2 

z=o . 

WV=W2 
4  W2=W V 

P  =  Y 

W1=DM-W2 
WW= 1 3 . 5*DM/W1 
R 1 =  98  7 . 7  +  4 . 992*WW 
R2=P**0 . 93022/77045.892 
C 

C  THE  PRESSURE  DROP  IN  THE  GIVEN  SECTION 

C  IS  EVALUATED  IN  SUBROUTINE  DPR. 

C  THE  RATE  OF  VAPOR  FORMATION  IS  EVALUATED 

C  IN  SUBROUTINE  DV. 

C 


. 


CALL  DPR (P, DP) 

CALL  DV ( PS , P , DP ) 

K 1 =H*DP 
L 1 =H*DW2 
Z  =  Z+H/2 . 

W2=WV+L1/2 . 

P=Y+K1/2 . 

R2=P**0. 93022/77045.892 

W1=DM-W2 

CALL  DPR ( P  f DP ) 

CALL  DV ( PS , P , DP ) 

K2=H*DP 
L2=H*DW2 
P=Y+K2/2 . 

W2=WV+L2/2 . 

R2=P** 0.93 02 2/7 7 04 5. 8 92 

W1=DM-W2 

CALL  DPR (P, DP) 

CALL  DV ( PS , P , DP ) 

K3=H*DP 
L3=H*DW2 
Z=  Z+H/2 . 

P=Y+K3 

W2=WV+L3 

R2=P**0. 93022/77045. 892 

W1=DM-W2 

CALL  DPR (P, DP) 


- 


. 


108 


CALL  DV ( PS , P , DP ) 

K4=H*DP 

L4=H*DW2 

Y=Y+ ( K 1 +2 . *K2+2 . *K3+K4 )/6. 

WV=WV+  (LI +2  .  *L2  +  2  .  *L3+L4 )/6. 

TEMP= 166. 766* Y* *0.06978-273. 

C 

C  TEST  TO  KNOW  WHEN  THE  TOP  OF  THE  SECTION 

C  IS  REACHED  (  Z*=  1  )  . 

C 

IF  ( Z . GE . 0 . 999 ) GO  TO  3 
GO  TO  4 
3  W2=WV 

W1=DM-W2 

C 

C  TEST  TO  KNOW  WHICH  SECTION  WAS  EVALUATED. 

C 

IF  ( NK . EQ . 2 ) GO  TO  5 
IF  ( NK . EQ . 3 ) GO  TO  6 
VAP=W2*66 . 

COND=W 1*66 . 

DEN=R2*U+ ( 1 . -U) *R1 
DP  1=0.  1 53*DEN*  9 . 8 1 
Y=Y-DP1 
GO  TO  7 

5  VAP=W2  *  111. 


COND=W 1  *  111. 


. 
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DEN=R2*U+ ( 1 ,-U) *R1 
DP2=0 . 1 53*DEN*9 . 8 1 
Y=Y-DP2 
GO  TO  7 

6  VAP=W2* 156. 

COND=W1*  1  56  . 

C 

C  OUTPUT  AT  THE  TOP  OF  EACH  SECTION: 

C  SECTION,  W1,  W2 ,  TEMP.,  ALPHA. 

C 

7  WRITE  ( 6 ,  1 0  0 ) NK , COND , VAP , TEMP , U 
100  FORMAT  (5X,I2,4( 1  OX, El  2. 5) ) 

IF  ( NK . LT . 3 ) GO  TO  8 
C 

C  CONVERGENCE  TEST  ON  THE  OUTLET  TEMPERATURE. 

C 

DT=88 . 5-TEMP 
TEST=ABS ( DT ) 

WRITE  (6,200) 

200  FORMAT  (/) 

IF  ( TEST . GT . 0 . 0 0 1 )  GO  TO  9 
STOP 


END 


' 


SUBROUTINE  DPR(P,DP) 


1  10 


c 

C  TO  EVALUATE  DP/DZ*  IN  ONE  TUBE  OF 

C  A  GIVEN  SECTION. 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , DW2 ,Y,Z,L,U,V1,V2 
REAL  J1 , J2 

IF  (W2.EQ.0. )  GO  TO  1 
A=3 . 141 6*R**2 
C= 1 . 20736E-05 
A 1  =  3  .  1  **2 

D=9.81*2.*R*(R1-R2) 

D=SQRT ( D ) 

C 

C  SOLVING  EQUATION  38  TO  DETERMINE  ALPHA. 

C 

J 1 =W 1 /SQRT ( R 1 ) /A/D 
J2=W2/SQRT ( R2 ) /A/D 
B 1 =  3 . 1*( J2+J1-1  .  ) 

C  1  =- J  1 

DELT=B 1 **2-4 . *A 1 *C 1 

W= ( -B 1 +SQRT (DELT ) ) /2 . /A  1 

U= 1 . -w 

S=W2*R 1 /W1/R2* ( 1 .-U)/U 
VI =W1/R1/A/( 1 . -U) 

V2=W2/R2/A/U 


C 


' 


1 1 1 


C  SOLVING  EQUATION  37. 

C 

RO=U*R2*S**2+ ( 1 . -U) *R1 
X=R2*U+R1*( 1 .-U) 

D=V 1 **2*RO*U*C/R2/P**0 . 06978 
D= 1 .-D 

B= (RO* (R1-S*R2 )/S/R1/R2+S- 1 . )*DW2*V1/A 
B=B+9 . 8 1 *L*X 
DP=-B/D 
GO  TO  2 

1  DP=- 9 . 8 1 *L*R 1 

2  CONTINUE 
RETURN 


END 


SUBROUTINE  DV(PS,P,DP) 
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C 

C  TO  EVALUATE  DW2/DZ* 

C 

COMMON  W 1 , W2 , R 1 , R2 , R , DW2 ,Y,Z,L,U,V1,V2 
REAL  LAM 

T= 1 66 . 766*P**0 . 06978 
LAM=3 1 35665 .7-2345. 714 *T 
CP  =  4 180  . 

A  =  3 . 1416*2. *R*166.766*L 
C 

C  DETERMINATION  OF  THE  HEAT  TRANSFER 

C  COEFFICIENT  (EQUATION  68). 

C 

IF  (U.GE. 0.619)  GO  TO  1 
COEFF=2539 . 

GO  TO  2 

1  COEFF  =  4657. 1-3423. 7*U 

2  CONTINUE 
C 

C  SOLVING  EQUATION  62. 

C 

DW2=A*COEFF* ( PS **0.06978-P* *0.06978) 

DW2=DW2- (W2/2 . +W 1 ) *CP* 1 66 . 766*0 . 06978/P** 0 . 93022*DP 

DW2=DW2 /LAM 

RETURN 


END 


. 

. 


